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Introduction 


Within the Crew Systems & Aviation Operations Branch (CSAOB) at LaRC, researchers are 
working to develop technology for use in the development of high-density air traffic operations 
for the evolving Next Generation Air Transportation System or NextGen as it is more commonly 
known. This technology is being developed within a project called Airspace Super Density 
Operations (ASDO) under NASA’s Airspace Systems Program (ASP). Researchers in CSAOB 
are working to develop real-time aircraft trajectory generation and real-time adjustment of the 
trajectories that results in energy-efficient trajectories for a given airplane model. A generated 
trajectory includes a lateral path profile and a vertical path profile where the range for the 
vertical profile is dependent on the lateral profile. To realize high-density operations, aircraft 
guidance and control systems need to be developed to accurately follow these real-time 
generated and real-time adjusted trajectories. As a result of this need, the authors — assigned to 
the Dynamic Systems and Control Branch (DSCB) at LaRC — developed an outer loop vertical 
navigation (VNAV) guidance and control system to track a vertical profile of a real-time 
generated trajectory. This document describes the VNAV guidance and control system, provides 
the specifications for the control system, presents VNAV guidance and control law gain 
determination, includes the stability analysis perfonned, and discusses the typical performance of 
the system. 

Overview of VNAV Outer Loop Control System 

The VNAV outer-loop control system for ASDO research was designed and developed to 
provide satisfactory control performance when coupled to a representation of the pitch inner loop 
control system that resides within a B757 flight control computer (FCC) as requested by 
CSAOB. For the B757 automatic flight control system architecture, the B757 VNAV outer loop 
control system resides in the flight management system (FMS) as illustrated in Figure 1 (note 
shaded block). The B757 VNAV outer loop control system outputs a pitch command signal 
(VNAVS) and a pitch rate command signal (VNAVI) for input to the FCC. The VNAV control 
laws described in this report were designed to provide compatible VNAVS and VNAVI signals 
to a representation of the B757 FCC that resides within a B757 aircraft simulation available to 
CSAOB. 

For the ASDO research, VNAV control laws were developed for path and speed control modes 
along with engagement logic. The control laws follow the ASDO trajectories using altitude and 
altitude rate errors, a speed command, control mode engagement requests (speed or path mode 
request), desired flight phase, and logic discretes provided by the ASDO VNAV trajectory 
calculations done in the Research Prototype Flight Management System (RPFMS) environment. 

The vertical trajectory calculations in the RPFMS compute a trajectory consisting of a series of 
connected straight line path segments and, in turn, compute altitude and altitude rate errors 
relative to the path segment that the aircraft is currently tracking or following (thus, referred to as 
the current path segment). To facilitate smooth transitions between segments, the trajectory 
calculations also compute altitude and altitude rate errors relative to the next path segment to be 
subsequently followed. These errors, along with altitude and altitude rate errors relative to a 
constraint altitude (an altitude that the aircraft should not fly through or away from is the selected 
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altitude shown on the mode control panel, or MCP, of the aircraft), are provided for the VNAV 
control system. 

The RPFMS trajectory calculations also compute a desired speed for each of the path segments 
of the VNAV trajectory. The desired speed is either a desired calibrated airspeed (CAS) or a 
desired Mach number (Mach). The trajectory calculations compute the desired true airspeed as a 
function of desired CAS or Mach number, and this airspeed is input to the VNAV control 
system. The computations in the RPFMS also send a discrete signal to the VNAV control 
system to indicate if the desired true airspeed input was computed from the desired CAS or the 
desired Mach number. 

When the VNAV control system is engaged in the path mode, it uses the RPFMS computed path 
errors for computing pitch steering signals to follow the VNAV trajectory and to level the 
aircraft flight at the constraint altitude should it be encountered. The VNAV control system 
outputs a pitch-steering signal (VNAVS) and a pitch-rate-steering signal (VNAVI), as shown in 
Figure 1, and these are input to the FCC for computing a pitch command (these signals will be 
described later in subsection VNAV Control System Outputs). The control law in the block 
Pitch Inner Loop uses the pitch command to compute the autopilot elevator servo command. 
When the VNAV control sytem is in the path mode, the pilot or an autothrottle system must 
maintain the desired CAS or desired Mach number by making appropriate throttle commands. 

When the RPFMS VNAV trajectory calculator requests the speed mode, the VNAV control 
system engages in the speed control mode where now the VNAVS and VNAVI output signals 
become signals for computing pitch commands in the FCC to track either desired CAS or desired 
Mach. In the speed control mode, the closed loop VNAV control system becomes a speed-on- 
elevator control mode; that is, the elevator changes the aircraft pitch to control the desired speed. 
When the VNAV control system is engaged in the speed control mode, the pilot or thrust 
management system is expected to control the throttle position for a fixed EPR — typically, either 
a climb EPR or an idle EPR setting depending on the phase of flight. 

The details within the blocks Pitch Command Process and Pitch Inner Loop of Figure 1 are 
proprietary to Boeing. However, the algorithm in Pitch Command Process combines and limits 
the VNAVS and VNAVI signals to produce a pitch command that will not cause a specified 
vertical acceleration limit to be exceeded. 
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Figure 1. Illustration of VNAV control laws location within a B757 flight control system architecture. 


Description of VNAV Control System 

Figure 2 illustrates the three basic parts that make up the VNAV control system: (1) 
Preprocessing for filtering sensor measurements and computing control law parameters, (2) 
Control Laws that compute the pitch and pitch rate steering output signals, and (3) Control Mode 
Logic for engagement and disengagement of the VNAV system control modes. These three 
parts and the inputs and outputs to the VNAV control system are described in the following 
subsections and defined in Tables 1 and 2. 
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Figure 2. VNAV outer-loop control system. 


VNAV Control System Inputs 

All of the VNAV control system inputs are defined in Table 1. The inputs include inertial 
vertical speed and flight path acceleration available from an aircraft’s inertial navigation system 
and true airspeed, Mach number, and barometric altitude available from an aircraft’s airdata 
computer. Some of these inputs are preprocessed or filtered to remove signal noise. The 
RPFMS trajectory calculations provide altitude and altitude rate error inputs relative to (1) the 
current VNAV path segment, (2) the next VNAV path segment, and (3) the constraint altitude. 
Other inputs that the RPFMS trajectory calculations provide are true airspeed command, 
requested control mode (path or speed), various arm/enable discretes, desired speed tracking 
discrete (MachSel), and current flight phase (climb, cruise, or descent). 
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Table 1. VNAV control system inputs. 


Parameter 

Units 

Description 

AIRCRAFT PARAMS 



Hdot 

ft/sec 

Inertial navigation system vertical velocity 

Mach 

Non-dim 

Mach number from air data computer 

VT 

ft/sec 

True airspeed from air data computer 

A fp 

ft/sec 2 

Inertial navigation system flight path acceleration 

NAV VARIABLES 



Ah, 

ft 

Altitude error relative to current path segment; h cmd i - h 

Alii 

ft/sec 

Vertical velocity error relative to current path segment; h cmcH - h 

Ah i+1 

ft 

Altitude error relative to next path segment 

Ah i+ i 

ft/sec 

Vertical velocity error relative to next path segment 

Ah c 

ft 

Altitude error relative to constraint altitude 

Ah c 

ft/sec 

Vertical velocity error relative to constraint altitude 

MachSel 

non-dim 

Mach select discrete: 0 - track commanded CAS, 1 - track 
commanded Mach number 

VSENABLE 

non-dim 

Vertical Speed (VS) submode discrete: 0 - inhibit vertical speed 
engagement, 1 - enable submode engagement 

DelTAS_VS 

ft/sec 

Constant true airspeed error value at which VS submode engages 

VNAV ARM 

non-dim 

VNAV mode arm discrete: 1 - armed, 0 - not armed 

PMODE 

non-dim 

Desired VNAV control mode: 0 - speed mode, 1 - path mode 

CurCapEnb 

non-dim 

Capture of current path segment enable discrete: 0 - disabled, 1 - 
enabled 

NSEGARM 

non-dim 

Next path segment arm discrete: 0 - not armed , 1 - armed (capture 
of next path segment armed) 

AltCapEnb 

non-dim 

Capture of constraint altitude enable discrete: 0 - disabled, 1 - 
enabled 

SEGUPD 

non-dim 

One-shot discrete set true for one iteration when aircraft arrives at 
start of a path segment 

FltPhase 

non-dim 

Phase of flight: 1 - climb, 2 - cruise, 3 - descent, 4 - altitude hold 

DISENG 

non-dim 

VNAV disengage discrete: 0 - do not disenage VNAV control, 1 - 
disengage VNAV control 

PARAMETERS 



PthHddLim 

ft/sec 2 

path control law error signal rate limit (vert accel limit) 

SpdHddLim 

ft/sec 2 

speed control law error signal rate limit (vert accel limit) 

VSHddLim 

ft/sec 2 

vertical speed control law error signal rate lim (vert accel limit) 

DHCP 

ft 

Path capture logic altitude error parameter 
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The altitude error inputs that are relative to path segments of the VNAV trajectory and provided 
by the RPFMS trajectory calculations are further clarified here. In this document the current 
path segment is referred to as the i th segment and the next path segment as the i th +l segment. 
Figure 3 depicts the altitude errors relative to the current and next path segments for two 
different scenarios. These two scenarios were chosen among several possible scenarios to 
illustrate that the path errors for the i th +l path segment are generated relative to the backward 
projection of the segment from its starting location. The path errors for the backward projection 
are used to trigger the capture of the path segment before arriving at the start of the segment so 
that the segment is captured without significant path overshoot. For scenario 1, the aircraft is 
below both the i th path segment and the backward projection of the i th +l segment (projected 
backwards since the starting point of the i th +l segment lies ahead of the aircraft); also, the flight 
path angle (y) of the i th +l segment is less than the flight path angle of the i th segment. For 
scenario 2, the aircraft is below the i th segment and above the backward projection of the i th +l 
segment; also, the flight path angle of the i th +l segment is greater than the flight path segment of 
the i th segment. The VNAV control system uses the altitude error and vertical velocity error 
relative to the i th path segment to compute pitch steering signals to track the i th path segment prior 
to capture of the i th +l segment. Prior to the capture of the i th +l segment, the VNAV control 
system uses the altitude error and vertical velocity error relative to the i th +l path segment to 
detennine when to capture this path segment so that a smooth and low-overshoot capture is 
achieved. At the instant when capture of the i th +l segment occurs, the errors relative to i th +l 
segment are then fed to the path control system for performing the capture maneuver. When the 
aircraft arrives at the intersection, the i th +l segment becomes the i th or current path segment. The 
capture logic is described later in the subsection Control Mode Engagement and Path Capture 
Logic. 
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Figure 3. Depiction of altitude errors for VNAV control system. 


VNAV Control System Outputs 

The control outputs from the VNAV control system are pitch and pitch rate steering commands 
that are suitable inputs for the pitch processing algorithms in the B757. The FCC processes these 
commands to produce a pitch command for the pitch inner loop control system. The FCC pitch 
processing algorithms compute a pitch error that is driven towards zero by the pitch inner loop 
autopilot elevator command control law. When the VNAV control system is tracking a path 
segment in path mode or the commanded speed in speed mode with zero error, the VNAV 
control steering commands will be zero. The other outputs of the VNAV control system are 
engagement discretes to indicate the status of the VNAV control system to the RPFMS trajectory 
computations. Table 2 provides a list of the VNAV control system outputs, their units, and 
descriptions. 
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Table 2. VNAV control system outputs. 


Parameter 

Units 

Description 

CONTROL COMMANDS 



VNAVI 

deg/sec 

VNAV pitch rate steering command 

VNAVS 

deg 

VNAV pitch steering command 

CONTROL DISCRETES 



VNAVE 

non-dim 

VNAV control mode engaged 

PTHARM 

non-dim 

VNAV path mode anned for engagement 

PTHENG 

non-dim 

VNAV path mode engaged 

SPDENG 

non-dim 

VNAV speed mode engaged 

VSENG 

non-dim 

VNAV vertical speed submode of speed mode 
engaged 

PTHCAP 

non-dim 

One-shot discrete true for one iteration when 
current path segment, next path segment, or 
constraint altitude is captured. 

TypePthCap 

integer 

0 - no path capture, 1 - current path segment 
capture, 2 - next path segment capture, 3 - 
constrained altitude capture 

NSEGCAP 

non-dim 

One-shot discrete true for one iteration when next 
segment path capture occurs; same infonnation 
available from above two outputs 


Preprocessing 

The preprocessing part of the VNAV control system consists of parameter filtering and control 
law parameter calculations. 

Parameter filtering is conducted to remove noise from the sensor measurements. The sensor 
measurements that are filtered are true airspeed (VT), flight path acceleration (A FP ), and Mach 

number (see Figure 4). The filtered airspeed (VT f ) is used in both of the VNAV control laws. 
True airspeed is further low-pass filtered to produce a “gain scheduled” filtered true airspeed 
(V T gs ) that is used in the calculation of control law parameters. Filtered flight path acceleration 

(A FPf ) and “gain scheduled” Mach number (Mach GS ) are both used to compute control law 
parameters. 

Figures 5 and 6 show that control law parameters are calculated for use in the VNAV mode 
control logic and VNAV speed control. The parameters are estimated weight-to-thrust change 
(EstDTW), climb discrete (CLMB), descent discrete (DESND), true airspeed error (AVT), and 
true airspeed rate-of-change bias (VdotB). The parameter AVT is used in the control mode logic, 
and the remaining parameters are used in the speed command processor (to be discussed in more 
detail in the Speed Control section of this document). EstDTW is used for computing 
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acceleration limits in the speed command processor, and its derivation is presented in Appendix 
A. The inputs used to compute EstDTW are the vertical velocity available from an inertial 
navigation system (Hdot) and two fdtered measurements discusssed above (VT GS and A FPf ). The 
true airspeed rate parameter, VdotB, is used in the speed command processor to reduce a 
potential bias offset in the processor output (Vcmd) and the derivation of this bias parameter is 
presented in Appendix B. More discussion of VdotB is given in the Speed Control section of 
this document. 
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Figure 5. VNAV parameter calculations. 



Figure 6. Algorithm for VNAV true airspeed rate-of-change parameter (VdotB). 


VNAV Control Laws 

Figure 7 shows the three control laws in the VNAV control system: 1) a VNAV path control law 
used for the path control mode, 2) a VNAV speed control law used for the speed control mode, 
and 3) a vertical speed control law. The vertical speed control law is a sub-mode of the speed 
control mode. All these control laws compute pitch and pitch rate steering commands that are 
output for processing by a pitch inner loop control system. 

When the path control mode is engaged, the VNAV path control law computes commands to 
track a desired VNAV trajectory, and speed is assumed to be controlled by the aircraft’s TMC 
autothrottle control law or by the pilot. When the speed control mode is engaged, one of the two 
control laws for this mode will be engaged. If the absolute value of the true airspeed error is less 
than or equal to a specified value (DelTASJVS of Table 1), the VNAV speed control law will 
track the desired VNAV speed profile. If the airplane is in the climb phase and the airspeed error 
becomes greater than DelTASVS (which occurs when the VNAV system commands a large 
airspeed increase or the pilot dials in a large airspeed increase on the MCP), then the vertical 
speed control law will engage and control to a constant vertical speed between 0 and 500 ft/min. 
When the true airspeed becomes less than DelTASJVS, the VNAV speed control law re- 
engages. Normally, the controlled vertical speed will be 500 ft/min since the airplane will 
normally be climbing at a higher rate than 500 ft/min when in the climb phase. The controlled 
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vertical speed can be a constant value between 0 and 500 ft/min if the commanded airspeed 
increase occurs while the airplane is climbing within that vertical speed range, or 0 if the airplane 
was temporarily descending when the airspeed increase was commanded. When climbing, the 
vertical speed control law insures that the airplane accelerates quickly to a large increase in the 
commanded airspeed, but without descending, as could result if the VNAV speed control law 
remained engaged. When descending, the vertical speed control law insures that the airplane 
decelerates quickly to a large decrease in the commanded airspeed without climbing, as could 
result if the VNAV speed control law remained engaged. When descending and the true airspeed 
error is less than -DelTAS VS, the vertical speed control law enages to control the vertical 
speed between 0 and -500 ft/min. When the airspeed error is greater than -DelTAS_VS, the 
VNAV speed control law re-engages. When engaged in the speed mode, the VNAV trajectory 
errors are ignored and the TMC or the pilot commands a fixed EPR or throttle setting for the 
climbing flight phase or an idle throttle setting for the descent flight phase. Figure 7 shows that 
the discretes PTHENG, SPDENG, and SPDE1 select the pitch steering and pitch rate steering 
commands that will be contained in the output signals VNAVS C and VNAVI C . 

Path Control Law 

The VNAV path control law was designed to track a VNAV trajectory that consists of several 
straight-line, connected path segments with various flight path angles. The control law captures 
and tracks a path segment with inputs of the altitude error, altitude rate error relative to that path 
segment, and fdtered true airspeed (see Figure 7). The control law architecture was adopted 
from the altitude capture and hold control mode of the Total Energy Control System (TECS) that 
was flight tested on a B737 research aircraft at Langley Research Center [1], The computed gain 
(KHERR) converts the altitude error (Ah) to a commanded change in vertical speed, and that 
commands an exponential reduction in the altitude error with a time constant equal to one 
divided by KHERR [2], The equation for the computation of KHERR and the limiting of the 
product of KHERR times Ah was taken from the TECS as implemented on the B737 research 
aircraft [3]. The rate of change of the vertical speed command is limited (see Figure 7) in order 
to limit the resulting vertical acceleration. The parameter PthHddLim (Table 1) specifies the rate 
limit of the command that is equivalent to a limit on vertical acceleration. The altitude rate error 
input (Ah) provides path damping. The division of the vertical velocity command and Ah by 
filtered true airspeed results in a flight path angle error relative to the air mass (if the vertical 
motion of the air mass is ignored) after the summing block in Figure 7. This error is multiplied 
by path control law pitch rate steering gain (KIVN) to produce a pitch rate steering command 
and mutiplied by the path control law pitch steering gain (KPVN) to produce a pitch steering 
command. KHERR is a function of PthHddLim, the altitude rate error, and the constant 
parameter KHBIAS. KHERR is also used in the path segment capture logic (to be discussed 
later) which results in proper capture initiation times for different values of PthHddLim. 

Engagement of this control law and capture of path segments is transient- free (i.e. no step 
command changes) since, at engagement and path capture, the pitch and pitch rate steering 
commands will both be zero. Control law engagement is discussed in the document subsection 
Engagement Logic, and the path capture logic is discussed in the section Path Capture Logic. 
The determination of the gains KIVN and KPVN will be discussed in the Control Law Gain 
Determination section of this document. 
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Speed Control 


The speed control law was designed to track a desired CAS or Mach number using pitch control 
of the aircraft (often referred to as speed-on-elevator control) and is in the middle level of Figure 
7. The input to the control law is computed true airspeed error (VTerr) that is the difference 
between the output of a second order speed command processor (Vcmd), to be discussed in detail 
in the following subsection, and filtered true airspeed. The input to the speed command 
processor is the selected true airspeed that the RPFMS computes from either the commanded 
Mach number or CAS, depending on whether commanded Mach or commanded CAS is to be 
tracked (note that the commanded Mach number or commanded CAS will be from either the 
VNAV trajectory commands or the speed selection on the MCP). In Figure 1 using selected true 
airspeed rather than Mach and CAS commands, only one control law is needed to track either 
CAS or Mach. The gain KVTERR controls the response time in reducing speed errors. The 
product of speed error and KVTERR is rate limited to a value specified by the parameter 
SpdHddLim to ultimately limit vertical acceleration during pitch changes for speed control. A 
washout filter on speed error with a 0.5 sec time constant is used to provide speed capture 
damping. The gain KCWO provides for adjustment of the damping. The sum of the rate limited 
speed error and damping signal is multiplied by gain KCIVN to produce a pitch rate steering 
command and multiplied by gain KCVN to produce a pitch steering command. This control law 
also provides transient-free engagement in that the pitch and pitch rate steering commands will 
be zero at engagement of the speed control law, since the true airspeed error is zero because the 
speed command proccessor output (Vcmd) is initialized to the true airspeed at control law 
engagement. The detennination of the gains KCIVN and KCVN will be discussed in the Control 
Law Gain Determination section of this document. 

The bottom left side of Figure 7 shows the vertical speed submode control law. The signal 
VSCMD is the vertical speed command that will have a value between 0 and 500/60 ft/sec when 
in the climb phase and between 0 and -500/60 ft/sec when in the descent phase. The first order 
filter at the left side of the control law serves two purposes. First, it is used to smoothly change 
the vertical velocity command from the instantaneous vertical velocity of the aircraft at 
engagement to the fixed commanded value. Second, the filter provides transient- free engagment 
of the control law since the vertical velocity error will be zero at engagement. The second first 
order filter is used to obtain acceptable gain and phase margins. The vertical velocity of the 
aircraft (Hdot) is subtracted from the output of the first filter to produce a vertical velocity error 
which in turn is divided by the filtered true airspeed and multiplied by 57.3 to compute a flight- 
path-angle-like error in deg. Multiplying this error by the pitch rate steering gain, KVSIVN, 
produces the pitch rate steering command. Multiplication of the output of the second filter by 
the pitch sterring gain, KVSVN, produces the pitch steering command. 

Figure 8 shows the computational algorithm for the vertical velocity command. The vertical 
velocity command depends on the flight phase inputs CLMB (climb phase) and DESND (descent 
phase) and the vertical velocity of the aircraft at the time the vertical velocity control law is 
engaged (VSENG). 
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Figure 8. Vertical velocity command for the vertical velocity control law. 


Second Order Speed Command Processor 

The speed command processor in Figure 9 computes a true airspeed error signal (VTerr) for 
input to the VNAV speed control law. The speed command processor smooths the selected true 
airspeed command (VT sel ) which is an input from the RPFMS trajectory calculator. As 
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mentioned earlier, VT sel can represent a commanded CAS or a commanded Mach depending on 

which is to be tracked. The use of the processor prevents large steps in the speed error from 
being introduced to the speed control law and limits the rate of change of the true airspeed 
command (Vcmd) that is output from the processor to rates that the aircraft can achieve. The 
processor consists of a second order, critically damped fdter that contains computed acceleration 
limits and speed limits. The bandwidth of the fdter was set to 0.12 rad/sec, resulting in a fdter 
output (Vcmd) of 80 percent of the value for a step input in approximately 25 sec and 98 percent 
in 50 sec, in the absence of any acceleration limits. This provides a reasonable time for a speed 
change to a new commanded speed. The rate of change of the fdter output is limited by 
computed acceleration limits. 

The acceleration limits are computed as a function of flight phase using the delta thrust to weight 
change computed parameter, EstDTW, as shown in Figure 9. When climbing, the upper 
acceleration limit is the product of EstDTW, acceleration due to gravity ( g ), and 0.6 and the 
lower limit is -1000/VT GS . As noted in the derivation for EstDTW in Appendix A, multiplying 

EstDTW by g provides a quantity consisting of the flight path acceleration plus an estimated 
potential flight path acceleration. Multiplying this product by 0.6 makes the upper acceleration 
limit conservative to what should be achieveable. The lower limit for the climb flight phase is 
not normally encountered since the aircraft is usually either maintaining a constant speed or 
accelerating. However, if a speed reduction during a climb is commanded, the lower limit 
quantity represents an achieveable deceleration level at the fdtered true airspeed, VT GS . When 

the airplane is descending, the lower limit is given as the product of -EstDTW, g, and 0.6 and the 
upper limit is 1000/VT GS . 

The Vcmd is the output of a limited integrator and the integrator is limited to values between 0 
and 1000 ft/sec. 



Figure 9. Speed command processor for VNAV control system. 
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The input VdotB shown at the lower middle portion of Figure 9 is used to offset a bias that 
would otherwise ocurr in Vcmd when the aircraft is climbing or descending under constant CAS 
or Mach. When the aircraft is being flown at constant CAS, true airspeed increases if the aircraft 
is climbing or decreases if descending. For constant Mach below 36,089 ft barometric altitude, 
true airspeed decreases if the aircraft is climbing and increases if descending. For constant Mach 
above 36,089 ft, true airspeed remains constant regardless of vertical speed since the atmospheric 
temperature is constant above this altitude. So, when true airspeed is changing for a constant 
CAS or Mach flight, the input to the second order fdter, VT sel will be changing, and a lag or bias 
would occur in Vcmd without the compensation that is provided by the input VdotB. Figure 6 
shows the computation of VdotB, and the derivation of it is in Appendix B. 

Control Mode Engagement and Path Capture Logic 

This section describes the logic to perform the following functions: 

• engage/disengage the path control, speed control, and vertical speed control laws 

• initialize the speed command processor 

• initialize three control law filters 

• initialize two control law rate limiters 

• detennine when to initiate the capture of a VNAV path segment 

• determine which set of altitude, and altitude rate errors to input to the path control law 

Engagement Logic 

Figure 10 shows the mode control logic that generates the discretes for engagement and 
initialization. The engagement logic determines which VNAV control mode should be engaged. 
The path mode is engaged and initialized (PTHENG is set true) if (1) the RPFMS is requesting 
that the pitch path control mode be engaged (Pth_R is true ), (2) the path mode is not engaged 
(the past value of PTHENG is false), (3) the control system is anned for engagement 
(VNAV ARM is true ), and (4) the discrete to capture a path segment has been set true (PTHCAP 
is true). When PTHENG is true, the path control law is engaged and its pitch and pitch rate 
steering commands are output from the VNAV control system. PTHENG is latched true and is 
reset false when either the speed mode is requested (Spd_R is true) or the disengage discrete 
(DISENG) is set true. The request for the speed or path mode to be engaged is computed from 
input variable PMODE (see Figure 11). PMODE is an input from the RPFMS trajectory 
calculations. If the speed control mode is requested (Spd_R is true) when the path control is 
engaged , the path control law will be disengaged, and, if VNAV ARM is true, the speed control 
mode will be engaged (SPDENG is set true) on the same computational iteration. If DISENG is 
set true, the VNAV control system will be disengaged. 

The bottom of Figure 10 shows that when the speed mode is engaged, the vertical speed mode 
can be engaged (i.e. VSENG can be set true). VSENG is set true for two conditions: (1) if the 
true airspeed error (AVT) is greater than DelTAS VS, the climb discrete (CLMB) is true, and the 
vertical speed mode is enabled for engagement (VSENB is true), and (2) if the airspeed error is 
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less than -DelTASVS, the descend discrete (DESND) is true, and VSENB is true. The discrete 
SPDE1 indicates the speed mode is engaged, but the vertical speed submode is not engaged. 
This discrete is used for selecting the correct output signals from the VNAV control system and 
for the initilization of two control law integrators, one fdter, and one rate limiter (Figure 7 and 
Figure 9). 


PTHENG p 



Figure 10. VNAV control mode engagement logic. 


Path Capture Logic 

Figure 1 1 shows the path capture logic (note that subscript P on some discretes represents the 
past value of that discrete). As previously stated, the VNAV path control law is designed to 
follow a vertical path that consists of a series of connected straight line path segments with 
various flight path angles (see the section titled Path Control Faw). The control law is designed 
to start the capture of a path segment with a transient-free control law engagement (that is, the 
pitch and pitch rate steering commands will be zero at engagement). The control law provides 
an exponential capture of a vertical path segment given the altitude error and altitude rate error 
relative to the path segment. The path capture logic determines which path segment should be 
captured and when a transition to the next segment should occur. If a current path segment 
capture should occur, CSEGCAP is set true; if a next path segment path capture should occur, 
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NSEGCAP is set true; and, if constraint altitude capture should occur, ASEGCAP is set true. 
These three discretes are set true for only one computational iteration by the One-Shot blocks. 
As shown for the computation of the variable TypePthCap in the lower right side of Figure 11, 
the capture discretes have a defined priority. The constraint altitude capture has the highest 
priority and the current path segment capture the lowest priority. For example, if current path 
segment and next path-segment captures were set true simultaneously, the next path segment 
capture would be executed rather than the current path segment capture. The discrete PTHCAP 
is set true when any type of capture occurs and is subsequently used in control law 
sychronization and output command smoothing during control mode changes or changing 
control to a different path segment. For any type of path capture to occur, the enabling discretes 
for each type of capture must be true (CurCapEnb for current path segment capture, NSEGARM 
for next path capture, and AltCapEnb for the constraint altitude capture). These enabling 
discretes are inputs from the RPFMS trajectory calculations. 

The CSEGCAP discrete will be set true if CurCapEnb is true, the past value of the path mode 
engage discrete (PTHENG p ) is false, the past value of the path mode engage request (Pth_R p ) is 

true, and the discrete output from the top block of Figure 1 1 labeled “Path Cap Trig” is true. 
Figure 1 1 shows the contents of this Path Cap Trig block. The contents of the other two Path 
Cap Trig blocks are the same. The computations for the gain KHERR in the Path Cap Trig block 
is the same as that for KHERR in the path control law (see Figure 7). Note that path capture will 
occur if the aircraft is within DHCP feet of the path segment regardless of what the altitude rate 
error is (the numerical value of this constant parameter is specified in Table 1 and is set in the 
RPFMS). DHCP is set to a small value to insure that path capture will occur if the aircraft is 
close to the desired path segment when the path mode is requested for engagement with the 
VNAV control system engaged in the speed mode or when the VNAV control system is anned 
for engagement and the path mode is requested to be engaged. 

The next path segment capture discrete (NSEGCAP) will be set true if the past value of the path 
mode engage request (Pth_R p ) is true, NSEGARM is true, VNAV ARM is true, and the discrete 

output from the middle Path Cap Trig block of Figure 1 1 is true. 

The constraint altitude capture discrete (ASEGCAP) will be set true if AltCapEnb is true, 
VNAV ARM is true, the past value of VNAV engage (VNAVE P ) is true, and the discrete output 
from the bottom Path Cap Trig block of Figure 1 1 is true. 


18 




When a path capture occurs, the appropriate set of altitude and altitude rate errors that are input 
to the VNAV control system from the RPFMS trajectory calculations must be selected for input 
to the path control law. Figure 12 shows the logic and switches for this selection. The discretes 
CNTLNSEG and CNTLASEG are used to set the switches for selection of the appropriate 
altitude and altitude rate error for the path control law. When CNTLNSEG is true, the control 
law is capturing the next path segment before the aircraft arrives at the starting point of this path 
segment (recall from Figure 3, that during initial capture of the next path segment, the control 
law is controlling to an extension of the next path segment). When the aircraft arrives at the 
starting point of the next path segment, the one-shot discrete path segment update (SEGUPD) is 
set true by the RPFMS trajectory calculations, and that causes CNTLNSEG to be reset to false. 
When SEGUPD is set true, the RPFMS trajectory calculations set the current altitude and 
altitude rate errors to those that previously were the next segment altitude and altitude rate errors 
since the aircraft has arrived at the starting point of the path segment that previously was the next 
path segment and now has become the current path segment. CNTLNSEG is also reset if 
ASEGCAP is set true. When CNTLASEG is set true, the control law is capturing or tracking the 
constraint altitude. CNTLASEG is reset to false if AltCapEnb is set false by the RPFMS and 
either NSEGCAP or CSEGCAP is set true or the past value of the speed mode request (Spd_R) 
is set true. 
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Figure 12. Path error selection logic. 


The control law output commands are smoothed to prevent potential steps in the commands 
when switching between control laws or when a path segment capture occurs. While the control 
law to be engaged will initially have zero commands at the switch time, the current control law 
may have non-zero commands because it may still be capturing a commanded target at switch 
time. When capturing a next path segment in the path mode, the initial commands for the next 
path segment will be zero, but the current path segment commands may be non-zero if the 
current segment capture has not been completed when next segment capture occurs. So, 
smoothing is provided by an easy-on/easy-off function. The smoothing and the logic to activate 
the smoothing is referred to as pitch command synchronization in this document. When pitch 
command synchronization is desired, the PSYNC discrete is set true to trigger the start of the 
synchronization. Figure 13 shows the logic for setting the PSYNC discrete. When PSYNC is set 
true, the previous pitch and pitch rate steering commands are stored, and then these stored 
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commands are linearly ramped out while the new VNAV pitch and pitch rate steering commands 
are linearly ramped in over a time period specified by the parameter FadeTime (default to 2 sec). 

The smoothing is done as follows: (1) let G be a gain that varies linearly from zero to one over 
the period FadeTime, (2) let VNAVSc and VNAVIc be the new pitch and pitch rate command 
signals, (3) let VNAVScp and VNAVIcp be the values of the stored past pitch and pitch rate 
command signals just before PSYNC becomes true, and (4) let VNAVS and VNAVI be the pitch 
and pitch rate command signals output from the VNAV control system. When PSYNC is set 
true, the output pitch and pitch rate command signals are computed as 

VNAVS = VNAVSc * G + VNAVScp * (1 - G) 

VNAVI = VNAVIc * G + VNAVIcp * (1 - G) 

When the value of G reaches one, G is reset to zero, and the output pitch and pitch rate command 
signals are computed as 

VNAVS = VNAVSc 

VNAVI = VNAVIc 



One-Shot 



Figure 13. Pitch command synchronization logic. 
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Control Law Gain Determination 


During initial control law development, the control law gains were chosen by examining the path 
and speed capture dynamic responses obtained from simulations using a non-linear B757 
Simulink’ simulation (Simulink is simulation software available from The MathWorks, Inc.) 
where the aircraft model part of the simulation is proprietary. A goal for the system design was 
to use fixed gains and avoid using gain scheduling. Such a goal seemed reasonable since the 
VNAV system is an outer-loop control system used during flight phases where the path and 
speed tracking response is desired to be relatively constant and slow. The system architecture is 
essentially commanding flight path angle changes from which a pitch command is computed. It 
was assumed that for slow changes in the commanded pitch, system stability would be 
maintained by the pitch inner loop gain scheduling as aircraft operating conditions changed. In 
attempting to use fixed gains, optimization algorithms available in the Simulink software 
environment were employed to find gains for the path and speed control laws that would 
minimize path and speed errors over the VNAV flight conditions. For optimization of the path 
control law gains, a cost function was used that had penalities for altitude error and pitch rate 
error. The penalty on pitch rate error avoided gains resulting in vertical path maneuvering 
uncomfortable for airplane passengers. The optimization algorithm [4, 5] detennined path 
control law gains while the VNAV control system was tracking a descending VNAV path 
trajectory. Optimal gains were determined at specified times along the trajectory at 5, 10, 15 sec, 
etc. The optimal gains found at the specified times did not necessarily have the same values. 
Fixed gains were computed by averaging the gains found at the specified times. Then, 
simulation runs were made with these fixed gains to see if the path capture and tracking 
performance was acceptable. Additionally, the Simulink Control Design Toolbox was used to 
determine the stability of the VNAV control system at 20 sec intervals along a descending 
VNAV trajectory using the fixed gains. For the speed control law, a similar process was used to 
determine optimal gains for a cost function having penalties on true airspeed error and pitch rate 
error. The same procedure as described above for the VNAV path control law was used to 
determine constant gains for the VNAV speed control law and their suitability. 

Table 3 shows the constant gains used for the VNAV control laws. The gains used for the 
VNAV speed control law are those determined by the optimization procedure. The gains used 
for the path control mode are approximately two times larger than those determined by the 
optimization procedure. Nichols plots for the path mode revealed that there was excessive gain 
margin with the optimized gains for the path mode. Since it was desired to improve the path 
tracking performance for some of the tested VNAV trajectories, the path control law gains were 
increased by a factor of approximately two from the optimized gains. Then, with these new 
gains, the gain and phase stability margins were checked and still found acceptable as is shown 
in the control system stability analysis section of this report. 
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Control Law Parameters & Gains 


Tables 3 and 4 list the gain values and control parameters for the VNAV control system. 
Table 3. Control law gains. 


Gain 

Units 

Value 

Description 

Path Control 

KIVN 

deg/sec/rad 

20 

Path control law pitch rate steering command gain 

KPVN 

deg/rad 

200 

Path control law pitch steering command gain 

Speed Control 

KVTERR 

non-dim 

1.0 

Speed control law gain on true airspeed error 

KCIVN 

deg/sec/ft/sec 

0.020014 

Speed control law pitch rate steering command gain 

KCVN 

deg/ft/sec 

0.13562 

Speed control law pitch steering command gain 

KCWO 

sec 

5.1992 

Speed control law washout fdter gain 

Vertical Speed Control 

KVSIVN 

sec 1 

0.3 

Vertical speed control law pitch rate steering command 
gain 

KVSVN 

non-dim 

1.0 

Vertical speed control law pitch steering command 
gain 


Table 4. Control law parameters. 


Parameter 

Units 

Value 

Description 

WV 

rad/sec 

0.12 

Bandwidth of speed control law seccond-order processor 

PthHddLim 

ft/sec 2 

1.6 

Vertical acceleration manuever limit for path control 

SpdHddLim 

ft/sec 2 

3.0 

Vertical accleration manuever limit for speed control 

VSHddLim 

ft/sec 2 

3.2 

Vertical acceleration manuever limit for vertical speed 
control 

DelTASVS 

ft/sec 

42.195 

magnitude of true airspeed error to engage vertical speed 
control 

FadeTime 

sec 

2.0 

Pitch command synchronization ramp-in/ramp-out time 

KHE 

1/sec 

0.08 

Max value of computed gain KHERR 

KHBIAS 

1/sec 

0.017 

bias value for computation of KHERR 


Control System Stability Analysis 

After the fixed gains were detennined as previously described, a stability analysis was perfonned 
to evaluate whether the VNAV control system gain and phase margins met the desired margins. 
The desired stability margins were a gain margin of at least 6 dB and a phase margin of at least 
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30 deg [6]. This analysis was accomplished using the Simulink Control Design Toolbox. These 
margins were computed at 20 sec intervals, while VNAV descending trajectories were simulated 
in the aircraft/VNAV guidance and control Simulink simulation. The descending trajectories 
started at 37,000 ft at Mach 0.7964 in level flight. At the end time of simulated descents, the 
aircraft altitude was around 6000 ft, and the aircraft calibrated airspeed was around 180 knots. 

Two 1100 sec simulated descents were made for the stability analysis. In one descent, the 
VNAV path was tracked by the control system engaged in the path tracking mode, and in the 
other descent, the VNAV speed profde was tracked by the control system engaged in the speed 
mode. For the path mode descent, the throttle was initially reduced to the idle position at a rate 
of -2 deg/s by engaging the autothrottle in the idle mode and at about 650 sec, the autothrottle 
was engaged in the speed mode to then track the VNAV trajectory commanded CAS. For the 
speed mode descent, the throttle was initially reduced to the idle position at a rate of -2 deg/s by 
engaging the autothrottle in the idle mode, and after that, the throttle remained at the idle position 
for the entire simulated time. 

As a 1100 sec descent simulation ran, the Simulink Control Design Toolbox stopped the 
simulation at specified times to compute the gain and phase of the closed-loop VNAV control 
system as a function of frequency. At each specified time, the toolbox was used to (1) compute a 
linear model about the current system operating condition (where the Tustin method was applied 
for computing the linear model), (2) compute the gain and phase as a function of frequency for 
the VNAV engaged control mode, (3) save the computed values, and then (4) continue the 
simulation to the next specified time. Nichols and Bode plots were then made from these saved 
gains and phases. Nichols charts were used to assess phase and gain margins for the VNAV 
control system. The Bode plots were used to detennine the system frequencies where the gain 
and phase margin might not meet the design requirements. These frequencies would then be 
used to design compensators to increase these margins to an acceptable level when gain tuning 
alone might not increase the margins to acceptable levels. 

The first specified time for both the VNAV path and speed control modes was at 10 sec, the next 
time at 21 sec, and then every 20 sec after that, unless this time occurred during the 2 sec time 
period after the discrete PSYNC was set (during this time period, pitch synchronization is 
occurring and these time periods were avoided to eliminate potential numerical problems in 
using the Simulink Control Design Toolbox). If a specified time occurred during the 2 sec pitch 
synchronization period, the linear model determination was done immediately after this time 
period. As a result of the specified stop times, 55 linear models for the VNAV path mode and 55 
linear models for the speed mode were computed for the closed-loop VNAV control system. 

Figure 14 shows the Nichols chart at 41 sec for the VNAV path mode. The linear model gain 
magnitude and phase were plotted for frequencies ranging from 0.01 to more than 39.8 rad/s. 
The diamond symbols on the -180 deg phase line mark the 6 dB gain margin points and the 
diamond symbols on the 0 db magnitude line mark the 30 deg phase margin points. The gain 
and phase margin requirements for the control system are satisfied with room to spare at this 
condition. 
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Nichols chart at t(second) = 41 

From: Gain/AfdsElevatorCommand To: FCC_pitch_inLp, (open loop) 



Phase (deg) 

Figure 14. Nichols chart for VNAV path mode linear model at 41 sec of simulated VNAV descent path. 


Figure 15 shows the Bode plot for time 41 sec which was made for a frequency range of 0.01 to 
100 rad/sec. From the Bode plot, the gain and phase for a specific system frequency can easily 
be read, which is not so from the Nichols plot. The gain and phase margins can also be 
determined from the Bode plot. For example, the gain margin is determined by the difference of 
the gain at -180 deg phase and 0 dB gain. As with the Nichols plot, this occurs in three places of 
the plot. Drawing a vertical line on the plot where the system phase is -180 deg helps to 
determine the gain margins. Likewise, the phase margin can be determined by the difference in 
the phase where the system gain is 0 dB and -180 deg phase. Notice for this plot that to 
determine where the system gain was unity (0 dB) would require magnifying the plot around the 
frequencies from 0.4 to 1 rad/s. However, for this stability analysis, the Nichols charts were used 
for determining the gain and phase margins while the Bode plots were made, as mentioned 
earlier, to show the system frequencies where the gain and phase margins might not meet the 
desired margins. 
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Bode plot: time = 41 seconds, VTOD 
KIVN = 20, KPVN = 200 



Figure 15. Bode plot for VNAV path mode linear model at 41 sec of VNAV simulated descent path. 


As mentioned earlier, 55 linear system models of the closed-loop system were computed for both 
the path control mode and speed control mode during the 1100 sec descent trajectory using the 
Simulink Control Design Toolbox. Gain and phase were computed to make Nichols charts and 
Bode plots for each of these linear models. Rather than including each individual plot in this 
report, one composite Nichols plot was made for all 55 linear models for the path control mode, 
and one composite Nichols plot was made for all 55 linear models of the speed control mode. 
Figure 16 shows the Nichols composite plot for the VNAV path control mode. As briefly stated 
at the start of this section, the aircraft conditions for these plots range from an altitude of 37,000 
ft to 6000 ft and the speed ranged from Mach 0.7964 to a CAS of 180 knots. As shown in the 
composite plot, the gain and phase margins of the control system path mode exceed the required 
gain and phase margins (indicated by the diamond symbols) by a large amount. 
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Figure 16. Composite Nichols plots for 55 operating conditions for path mode of VNAV control system. 

Figure 17 shows the composite Nichols plot of gain and phase for the 55 linear models for the 
speed mode of the VNAV control system. These plots show that the gain and phase margins for 
all of the 55 linear models exceed the required gain and phase margins indicated by the diamond 
symbols on the composite plot. For some linear models, the gain margin was around 10 dB in 
the low frequency range as shown by the traces in the vicinity of the top diamond symbol. The 
conditions for which the gain margin was around 10 dB occurred between 461 and 601 sec 
where the altitude ranged from approximately 19,000 ft to 10,000 ft, and CAS was constant at 
300 knots. To detennine the approximate system frequency for this gain margin, a Nichols plot 
was made for one of those conditions at 501 sec and is shown in Figure 18. The plot shows that 
the gain margin of approximately 10 dB above the top diamond symbol occurs for a system 
frequency a little less than 0.16 rad/s, which is more than twice the phugoid frequency of the 
aircraft model at this speed. 
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Figure 17. Composite Nichols plots for 55 operating conditions for speed mode of VNAV control system. 
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Nichols chart att(second) = 501 

From: Gain/AfdsElevatorCommand To: FCC_pitch_inLp, (open loop) 



Phase (deg) 

Figure 18. Nichols plot of linear model at 501 sec for VNAV control system speed mode. 


Control Performance 

As mentioned previously, the control system was developed using a six DOF, non-linear B757 
Simulink simulation. The performance of the control system shown in this section is from data 
generated with that simulation. Trajectory data was provided by CSAOB in the form of 
spreadsheet files and are trajectories generated by CSAOB trajectory generation software. 
Simulink code was written to use that data to generate the trajectory errors and to emulate the 
RPFMS outputs needed by the VNAV control system. As mentioned in the VNAV Control 
Laws section of this report, an autothrottle system was used to control the speed of the aircraft 
when engaged in the VNAV path control mode. In order to determine the performance of the 
path control mode with the Simulink simulation, an autothrottle, based on an autothrottle control 
law used on the research B737 aircraft [1], was implemented in the simulation. Appendix C 
describes this autothrottle. 

The VNAV trajectory data defines the trajectory as connected straight path segments. The data 
provided for each path segment are the tangent of the flight path angle of the segment, ground 
range from a specified origin to starting point of the segment, altitude at the starting point of 
segment, flight phase (climb - 1, cruise - 2, descent - 3), desired pitch control mode (path mode 
- 1 or speed mode - 0), throttle control mode (fixed - 1, speed control - 2, or idle - 3), desired 
CAS in knots, desired Mach number, desired flap handle position in degrees, and desired gear 
position (up - 0, down - 1). Spreadsheet data files were generated from the trajectory data files 
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provided by CSAOB and a trajectory data fde was generated to specifically test captures of path 
segments with various flight path angles and to test switching between path and speed control 
modes. The two data files generated from the CSAOB data were named vnav data LAX- 
ATL_TOD.xls and vnav_data_LAX-ATL_CLMB.xls. The data file to specifically test path 
segment captures was named vnavCLMBdatal.xls. Tables 5-7 show these data files with text 
added to the top of the columns to identify the data. These data files were used to determine the 
performance of the VNAV control system for the path control and speed control modes. 


Table 5. Trajectory data file vnav_data_LAX-ATL_TOD.xls. 


Seg num 

Tangent fit 
path angle 

Ground 
range, ft 

Altitude, ft 

Fit phase 

Pitch m ode 

Thrust m ode 

CAS, knots 

Mach 
num ber 

Flaps 

Gear pos 

1 

0 

0 

37000 

2 

1 

2 

300 

0,7964 

0 

0 

2 

-0.0614876 

16000 

37000 

3 

1 

3 

300 

0,7964 

0 

0 

3 

-0.0678675 

26329,3963 

36364,87 

3 

1 

3 

300 

0,7964 

0 

0 

4 

-0.0734726 

55798,5565 

34364,87 

3 

1 

3 

300 

0,7964 

0 

0 

5 

-0.0781847 

83019,554 

32364,87 

3 

1 

3 

300 

0,7964 

0 

0 

6 

-0.0572472 

108600 

30364,87 

3 

1 

3 

300 

0,7964 

0 

0 

7 

-0.0573330 

141107,218 

28503,92 

3 

1 

3 

300 

0,7964 

0 

0 

8 

-0.0577469 

162738,189 

27263,75 

3 

1 

3 

300 

0,7964 

0 

0 

9 

-0.0582579 

197372,048 

25263,75 

3 

1 

3 

300 

0,7964 

0 

0 

10 

-0,0592010 

231702,1 

23263,75 

3 

1 

3 

300 

0,7964 

0 

0 

11 

-0,0586302 

265485,302 

21263,75 

3 

1 

3 

300 

0,7964 

0 

0 

12 

-0,0607301 

268219,554 

21103,44 

3 

1 

3 

300 

0,7964 

0 

0 

13 

-0,0619882 

301152,1 

19103,44 

3 

1 

3 

300 

0,7964 

0 

0 

14 

-0,0623943 

333416,274 

17103,44 

3 

1 

3 

300 

0,7964 

0 

0 

15 

-0.0631780 

345021,654 

16379,33 

3 

1 

3 

300 

0,7964 

0 

0 

16 

-0 .0640384 

376678,216 

14379,33 

3 

1 

3 

300 

0,7964 

0 

0 

17 

-0,0641701 

407909,45 

12379,33 

3 

1 

3 

300 

0,7964 

0 

0 

18 

-0,0159142 

418178.085 

11720,39 

3 

1 

3 

230 

0,7964 

0 

0 

19 

-0,0399531 

463445.145 

11000 

3 

1 

3 

230 

0,7964 

0 

0 

20 

-0,0398647 

535264,831 

8130,58 

3 

1 

3 

230 

0,7964 

0 

0 

21 

-0,0399894 

557624,936 

7239,2 

3 

1 

3 

210 

0,7964 

0 

0 

22 

-0,0231346 

588613,125 

6000 

3 

1 

2 

210 

0,7964 

0 

0 

23 

-0,0226698 

597909,581 

5784,93 

3 

1 

2 

206 

0,7964 

0 

0 

24 

-0,0231653 

600765,356 

5720,19 

3 

1 

2 

186 

0,7964 

1 

0 

25 

-0,0234017 

611702,364 

5466,83 

3 

1 

2 

180 

0,7964 

5 

0 

26 

-0,0590149 

614558,138 

5400 

3 

1 

2 

180 

0,7964 

5 

0 

27 

-0,0590574 

624583,728 

4808,34 

3 

1 

2 

177,2 

0,7964 

5 

0 

28 

-0,0525868 

641657,613 

3800 

3 

1 

2 

166 

0,7964 

5 

0 

29 

-0,0530980 

660918,899 

2787,11 

3 

1 

2 

163,4 

0,7964 

20 

0 

30 

-0,0532455 

662559,45 

2700 

3 

1 

2 

146 

0,7964 

20 

0 

31 

-0,0523814 

673071,13 

2140,3 

3 

1 

2 

138 

0,7964 

25 

1 

32 

-0,0532617 

674833,204 

2048 

3 

1 

2 

138 

0,7964 

25 

1 


30 



Table 6. Trajectory data file vnav_data_LAX-ATL_CLMB.xls. 


Seg num 

Tangent fit 
path angle 

Ground 
range, ft 

Altitude, ft 

Fit phase 

Pitch m ode 

Thrust m ode 

CAS, knots 

Mach num ber 

Flaps 

Gear pos 

1 

0.0222008 

0 

2000 

1 

0 

1 

210 

0.7964 

0 

0 

2 

0.0211934 

1154.46195 

2025.63 

1 

0 

1 

218 

0.7964 

0 

0 

3 

0.1489184 

4374.80316 

2093.88 

1 

0 

1 

240 

0.7964 

0 

0 

4 

0.1430454 

14218.1103 

3559.73 

1 

0 

1 

240 

0.7964 

0 

0 

5 

0.1360156 

25215.8793 

5132.91 

1 

0 

1 

240 

0.7964 

0 

0 

6 

0.1270879 

39920.0788 

7132.91 

1 

0 

1 

240 

0.7964 

0 

0 

7 

0.1189206 

55657.218 

9132.91 

1 

0 

1 

240 

0.7964 

0 

0 

8 

0 

62948.5566 

10000 

1 

0 

1 

240 

0.7964 

0 

0 

9 

0 

66898.0317 

10000 

1 

0 

1 

254.8 

0.7964 

0 

0 

10 

0,0896887 

81784.5146 

10000 

1 

0 

1 

308 

0.7964 

0 

0 

11 

0.0834758 

104083.859 

12000 

1 

0 

1 

308 

0.7964 

0 

0 

12 

0.0810733 

111375.197 

12608.65 

1 

0 

1 

308 

0.7964 

0 

0 

13 

0.0748086 

136044.226 

14608.65 

1 

0 

1 

308 

0.7964 

0 

0 

14 

0.0676688 

162779.134 

16608.65 

1 

0 

1 

308 

0,7964 

0 

0 

15 

0.0658315 

171711.024 

17213.06 

1 

0 

1 

308 

0.7964 

0 

0 

16 

0.0588317 

202091.602 

19213.06 

1 

0 

1 

308 

0.7964 

0 

0 

17 

0.053872 

227915.092 

20732,3 

1 

0 

1 

308 

0.7964 

0 

0 

18 

0.0473608 

265040.158 

22732.3 

1 

0 

1 

308 

0.7964 

0 

0 

19 

0.0411962 

307269.161 

24732.3 

1 

0 

1 

308 

0.7964 

0 

0 

20 

0.0355078 

355817.324 

26732.3 

1 

0 

1 

308 

0.7964 

0 

0 

21 

0.0313648 

412142.914 

28732.3 

1 

0 

1 

308 

0.7964 

0 

0 

22 

0.0433673 

425328.085 

29145.85 

1 

0 

1 

308 

0.7964 

0 

0 

23 

0.0407373 

471445.802 

31145.85 

1 

0 

1 

308 

0.7964 

0 

0 

24 

0.0373195 

520540.815 

33145.85 

1 

0 

1 

308 

0.7964 

0 

0 

25 

0.0307305 

574132.154 

35145.85 

1 

0 

1 

308 

0.7964 

0 

0 

26 

0 

634467.981 

37000 

2 

1 

2 

308 

0.7964 

0 

0 

27 

0 

890393.965 

37000 

2 

1 

2 

308 

0.7964 

0 

0 
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Table 7. VNAV trajectory file vnavCLMBdatal.xls. 


Seg num 

Tangent fit 
path angle 

Ground 
range, ft 

Altitude, ft 

Fit phase 

Pitch m ode 

Thrust m ode 

CAS, knots 

Mach num ber 

Flaps 

Gear pos 

1 

0.0524078 

0 

5000 

1 

1 

2 

205 

0.7964 

1 

0 

2 

0 

4961.1 

5260 

1 

1 

2 

205 

0.7964 

1 

0 

3 

0.185339 

40000 

7000 

1 

1 

2 

200 

0.7964 

1 

0 

4 

0.0524078 

45395.5 

8000 

1 

1 

2 

205 

0.7964 

1 

0 

5 

0 

90000 

5000 

1 

1 

2 

240 

0.7964 

0 

0 

6 

0.0524078 

96000 

5000 

1 

1 

2 

240 

0.7964 

0 

0 

7 

-0.0524077 

130000 

6000 

1 

1 

2 

240 

0.7964 

0 

0 

8 

0 

136000 

5680 

1 

1 

2 

240 

0.7964 

0 

0 

9 

0 

160000 

37000 

1 

1 

2 

300 

0.7974 

0 

0 

10 

-0.0699268 

171000 

37000 

1 

1 

2 

300 

0.7974 

0 

0 

11 

0 

230000 

37000 

1 

1 

2 

300 

0.7974 

0 

0 

12 

-0.0524777 

241000 

37000 

1 

1 

2 

300 

0.7974 

0 

0 

13 

0.0561258 

300000 

11000 

1 

1 

2 

250 

0.7964 

0 

0 

14 

0.0854929 

310000 

11561.3 

1 

0 

1 

311 

0.7974 

0 

0 

15 

0 

400000 

2000 

1 

1 

2 

210 

0.7974 

0 

0 

16 

0.0524078 

405000 

2000 

1 

0 

1 

250 

0.7964 

0 

0 

17 

0.0104724 

450000 

11000 

1 

1 

2 

250 

0.7964 

0 

0 

18 

0.0524078 

455000 

11000 

1 

0 

1 

311 

0.7964 

0 

0 

19 

0 

500000 

37000 

2 

1 

2 

311 

0.7964 

0 

0 


Path Control Performance 

Four Simulink simulation runs are included in this report to show the path capture and tracking 
performance of the VNAV control system. The runs consisted of (1) a top-of-descent (TOD) 
260 sec run from a cruise altitude of 37,000 ft, (2) an idle descent 400 sec run from 16,360 ft, (3) 
a 400 sec descent run from 7,235 ft to 1,260 ft with speed changes that resulted in flap 
deployment to 20 deg, and (4) a 200 sec run to show cruise altitude capture performance. For 
each run, the report includes at least three figures that are plots of various variables and logic 
discretes. Figures 19-21 show the performance for the first simulation run and will be used to 
explain the plotted variables and discretes. Figure 19 shows plots of altitude, altitude rate 
(Hdot), pitch attitude, pitch rate, vertical acceleration, flight path acceleration, elevator position, 
stabilizer position, commanded CAS, CAS, commanded power lever angle (PLA) or commanded 
throttle, and PLA. Figure 20 shows plots of altitude error (altitude command minus altitude) and 
altitude rate error (altitude rate command minus altitude rate) relative to the VNAV trajectory 
path segments, the pitch steering command (VNAVS), and the pitch rate steering command 
(VNAVI) output from the VNAV control system, the pitch attitude error (DTHETE or pitch 
attitude minus pitch command), the elevator servo command (DESCMD), the CAS error, and the 
Mach error. Figure 21 consists of six plots showing various VNAV discretes. The discretes are 
separated on the plots by an integer multiple of two and the false value of a discrete is indicated 
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by even numbers. True values of the discretes are an increment of one above the false values. 
The discretes for each plot are plotted in the order shown in the plot legends where the discrete at 
the top of the legend is plotted at the top of the plot. The first plot at the top of Figure 21 shows 
four discretes: the discrete to ann the VNAV control system (VNAV ARM), the discrete that 
requests engagement of the path control mode (VNAV_Pth_R), the discrete that requests 
engagement of the speed control mode (VNAV_Spd_R), and the discrete to disengage the 
VNAV control system (VNAV_Diseng). The second plot shows five discretes: the discrete that 
shows when the VNAV control system was engaged (VNAVE), the discrete that shows when the 
path mode was engaged (VNAV_Pth_Eng), the discrete that shows when the speed mode was 
engaged (VNAV_Spd_Eng), the discrete that shows when the vertical speed submode of the 
speed mode was engaged (VNAV_VS_Eng), and the discrete that shows if Mach number or 
CAS was being tracked (MACHSEL). The third plot shows six discretes related to path capture: 
the one-shot discrete that shows when smoothing of the pitch steering output commands started, 
the one-shot discrete that shows when capture of the current path segment occured, the one-shot 
discrete that shows when capture of the next path segment occured, the one-shot discrete that 
shows when capture of the constraint altitude occured, the one-shot discrete that indicates when a 
path segment capture occured (PTHCAP), and the discrete that showed when the control system 
was armed for capture of the next path segment (NSEG ARM). The fourth plot shows three 
discretes that show the requested thrust or throttle control mode: the discrete that shows when the 
constant thrust mode (THR R) was requested, the discrete that shows when the airspeed control 
mode (SPD R) was requested, and the discrete that shows when idle throttle was requested 
(IDLE R). The fifth plot shows VNAV control mode requests from the RPFMS system: the 
discrete that shows when the RPFMS requested the path control mode to be engaged 
(Pth mode R), and the discrete that shows when the RPFMS requested the speed control mode 
to be engaged (SpdmodeR). The sixth plot shows three discretes: the discrete that shows when 
the VNAV control system was capturing or tracking the next VNAV path segment but has not 
reached the starting point of this path segment (CNTLSEG), the discrete that shows when the 
VNAV control system was capturing or tracking the constraint altitude (CNTLASEG), and the 
one-shot discrete from the RPFMS trajectory calculations that shows when the aircraft has 
arrived at the starting point of a VNAV path segment (segupd) and that the RPFMS trajectory 
calculations has switched from what was the next segment trajectory errors to now be the current 
path segment trajectory errors. 

The first simulation run to show the path control performance is a TOD 260 sec run. The run 
begins with the aircraft flying level at 37,000 ft at Mach number 0.7964. The VNAV control 
system proceeds to capture a descending path segment and, at the same time, the engine throttles 
(PLA position) are open-loop driven to the idle position (PLA = 55 deg) at a rate of -2 deg/s. 
The initial aircraft weight for this run was 172,841 lb with a center of gravity (CG) located at 25 
percent of the mean aerodynamic chord (MAC). Table 5 shows the VNAV trajectory data file 
for this run. As mentioned above, fig. 19-20 show the plots of variables and discretes for this 
run. At 0.22 sec the VNAV control system is anned, the VNAV control system is engaged in the 
path mode, and both a current path segment capture (path segment 1) and next path segment 
(path segment 2) occurr. The one-shot discretes for current path capture (CSEGCAP) and next 
path segment capture (NSEGCAP) cannot be seen on the third plot of Figure 21 because of the 
time scale (the discretes are hidden by the y-axis). Given that next path segment capture has 
priority, the path control system begins to pitch the aircraft down to capture and subsequently 
track path segment 2. For all the path control runs, the vertical acceleration parameter 
PthHddLim was set to 1.6 ft/sec 2 (or 0.05 g). The fifth plot of Figure 19 depicts the vertical 
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acceleration which shows that a couple of peak values in the time period from 0 to 25 seconds 
exceeded 1.6 ft/s 2 - some with values near 2 ft/s 2 . The aircraft pitch response is shown to be 
smooth, as would be expected from the small values shown for pitch rate. Figure 20 shows the 
altitude and altitude rate error (Hdot error) histories. The altitude error jumps when a next 
segment capture occurs and some capture overshoot occurred. After path segment capture was 
completed (when Hdot becomes relatively constant), altitude errors are generally less than 5 ft 
and altitude rate errors are less than 1 ft/s. Since the throttle was driven to idle for this run and 
the control system was tracking the path, desired Mach and CAS were not actively followed. 
The VNAV trajectory was calculated so that for the first 140 sec of the run, the aircraft speed 
would follow a constant Mach and, after that time, the aircraft would follow a constant CAS (see 
Mach select discrete MachSel on Figure 21). The run ended with the VNAV control system 
tracking path segment 9 of the trajectory in Table 5. This path segment was captured at 230 sec. 
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Figure 19. Various longitudinal parameters for TOD run from 37,000 ft in path control mode. 
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Figure 20. Path errors and command signals for TOD run from 37,000 ft in path control mode. 
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Figure 21. Logic discretes for TOD run from 37,000 ft in path control mode. 
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The second simulation run showing VNAV control system path control performance is a 400 sec 
idle descent run from 16,360 ft. The aircraft is initially trimmed with a CAS of 300 knots at a 
flight path angle of -3.6 deg which is slightly greater than the flight path angle of path segment 
15 of the VNAV trajectory in Table 5 (-3.615 deg). The aircraft trim weight is 172,690 lb with 
the CG located at 19 percent MAC. Figures 22-24 show time histories for this run. The current 
path segment (path segment 15) is captured around 10 sec as shown by the one-shot discrete 
CSEGCAP in Figure 24 and, at that time, the VNAV control system is engaged in path mode as 
shown by discrete VNAV_Pth_Eng in Figure 24, and the throttle begins to reduce to idle at -2 
deg/s. The capture of path segment 16 occurs around 35 sec. At 100 sec, path segment 18 is 
captured and it has a signficantly shallower flight path angle than path segment 17, as is 
illustrated by the vertical velocity (Hdot) in Figure 22. Figure 22 shows that the shallower flight 
path angle results in CAS decreasing from 300 knots towards the desired speed of 230 knots. 
The maximum vertical acceleration in capturing path segment 18 is 2.1 ft/s 2 that exceeded the 
desired maximum of 1.6 ft/s 2 specified by parameter PthHddFim. Around 180 sec, a next 
segment capture of path segment 19 occurs which has a steeper flight path angle than segment 18 
to maintain CAS at 230 knots. However, Figure 22 shows that CAS continues to decrease below 
230 knots at 230 sec. The last path segment captured is segment 21 at around 370 sec, and 
shortly after that, the flaps are command to 1 deg. The deployment of the flaps to 1 deg causes a 
small wiggle in the pitch rate. The plot of altitude error in Figure 23 shows there is a 10 ft 
overshoot in capturing path segment 18 and about a 5 ft undershoot in capturing path segment 
19. When tracking the path segments the altitude error is less than 10 ft. The small wiggles in 
the pitch rate in Figure 22 from 150 sec to 340 sec are due to trimming of the stabilizer position. 
Trimming occurs when the elevator control deflection persists at some deflection away from the 
elevator neutral position (the neutral position is the position that results is zero force at the pilot’s 
control column). Trimming stops when the elevator is returned the the neutral position. When 
the stabilizer is trimming the elevator is driven in an opposite direction in an attempt to cancel 
the pitch moment generated by the changing stabilizer position. The amount of elevator 
deflection to cancel pitch moment generated by the stablizer movement has to be estimated from 
aircraft measurements using table lookup data. This estimate is not perfect resulting in the small 
wiggles in pitch rate. Figure 22 shows the stabilizer movement along with the elevator surface 
position. 
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Figure 22. Various longitudinal parameters for idle descent from 16,360 ft in path control mode. 
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Figure 23. Path errors and command signals for idle descent from 16,360 ft in path control mode. 
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Figure 24. Logic discretes for idle descent from 16,360 ft in path control mode. 
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Figures 25-28 show the results of the third 400 sec simulation run to show path control 
performance of the VNAV control system. Figure 25 includes a fourth plot to show the various 
changes in the flap position that occur in this run along with some variables. The simulation run 
begins with the aircraft trimmed at an altitude of 7,235 ft with a flight path angle of -2.2 deg, 
CAS of 210 knots, flaps up, and the aircraft flying toward path segment 21 of the VNAV 
trajectory in Table 5. The VNAV control system engages immediately when it is armed for 
engagement (VNAV ARM set true) at 0.22 sec. At this time, the control system begins capture 
of the current path segment. The discrete CSEGCAP cannot be seen in Figure 27 because it is 
hidden in the y-axis because of the time scale for this run. The requested idle thrust mode for 
path segment 21 is over-ridden for this run and set to the speed control mode (Figure 27). The 
throttle (PLA position) remains close to the trim value of 60 deg until path capture of segment 24 
at 103 sec when CAS is commanded to 186 knots. Then, the throttle or PLA moves close to the 
idle value of 55 deg. When tracking path segment 21, the altitude error is essentially zero and 
CAS error is less than 1 knot (Figure 26). The capture of path segment 22 occurs around 67 sec 
as a next-path- segment capture, and, at around 79 sec, while the capture of path segment 22 is 
underway, the capture of path segment 23 ocurrs. The relatively small time between path 
captures occurs because path segment 22 is very short (approx. 9,300 ft in ground range). 
Around 140 sec, path segment 26, whose flight path angle is steeper than that of path segment 
25, is captured. The vertical acceleration reaches about -2.1 ft/s 2 which exceeds the desired 
maximum magnitude of 1.6 ft/s 2 . When path segment 24 is captured and CAS is commanded to 
186 knots around 105 sec, the flaps are commanded to 1 deg (Figure 28). Although the desired 
flaps for the path segment is 1 deg (Table 5), the flaps are not deployed using the desired flap 
position from the trajectory data. Instead, the flaps are commanded from a Simulink subsystem 
block to various positions according to a schedule as a function of CAS reference speeds and 
aircraft weight for flaps 25 and 30. As the various path segment captures occur, CAS is 
commanded to a smaller value (Figure 25 and Figure 28). The commanded CAS is not tracked 
closely after capture of path segment 28 around 220 sec, because the PLA is essentially at idle 
position. However, CAS continues to decrease for this run, and is within 5 knots of the 
commanded 138 knots for path segment 31 and is the last path segment captured at 306 sec for 
this simulation run. Ligure 28 shows that the flaps are commanded to 15 at 350 sec and to 20 at 
370 sec. Ligure 26 shows that the flap deployment causes the aircraft to rise above the desired 
path by 20 ft at 375 sec (negative altitude error is aircraft above the desired path). However, the 
control system is then correcting back to the path having an altitude error of -4 feet at the end of 
the run and still converging. Note that the vertical acceleration reaches 2 ft/s 2 which exceeds the 
desired maximum magnitude of 1.6 ft/s 2 . In general, most of the time, the altitude error after 
path segment capture while tracking the path segments is less than 5 ft. The aircraft altitude at 
the end of the run was 1,261 ft which is less than the 2000 ft that is specified for the constraint 
altitude. The aircraft flys through the constraint altitude because constraint altitude capture is not 
enabled for this run. 
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Figure 25. Various longitudinal parameters for descent from 7,235 ft in path control mode. 
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Figure 26. Path errors and command signals for descent from 7,235 ft in path control mode. 
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Figure 27. Logic discretes for descent from 7,235 ft in path control mode. 


45 






40 


a: 30 

Qj 

T> 

a 20 

as 

cr io 


1 1 1 1 I — 

Flap Hndl 

flaps 


/ 




i 


r 



Figure 28. Longitudinal variables for descent from 7,235 ft in path control mode. 
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Figures 29-30 show the results for the fourth simulation run showing cruise altitude capture 
performance. The run begins with the aircraft trimmed at an altitude of 35,146 ft with a flight 
path angle of 1.7613 deg, Mach number of 0.7964, weight of 195,906 lb, and the CG located at 
20 percent MAC. The initial position of the aircraft is at the start of path segment 25 for the 
trajectory in Table 6. The VNAV control system is armed for engagement at 0.22 sec and 
engages at that time in the speed control mode as that control mode is being requested (see 
Figure 31 - note VNAV Spd Eng and SpdmodeR discretes). The control system tracks the 
Mach number since the discrete MACHSEL is true (Figure 31). At 20 sec, the path mode is 
requested to be engaged (see VNAV Pth R and Pth mode R discretes in Figure 31). The path 
mode does not immediately engage and capture path segment 25 because the current path capture 
enable discrete is set false for this run. The next level cruise altitude segment 26 is not captured 
when the path mode is requested because the aircraft is not close enough for capture at this time. 
Thus, the control system remains in the speed control mode until 69 sec when the next path 
segment capture discrete is set true (see NSEGCAP in Figure 31). Then, the control system 
disengages from the speed mode and engages in the path mode to begin capture of the level 
cruise path segment 26. The capture of the cruise altitude is overshot by 4 ft, and, after that, the 
altitude error is reduced. However, the plot of altitude error in Figure 30 shows that there is an 
approximately 20 sec period low-amplitude stable oscillation in the altitude error with an 
amplitude between ±1 ft after 165 sec. Lowering the control gains did not eliminate this 
oscillation. The oscillation was determined to be caused by the 0.06 deg hysteresis on the 
elevator command that is contained in the elevator control surface actuator model. This 
hysteresis is evident in the plot of the average left and right elevator surface positions of Figure 
29. This plot shows after 110 sec that the average elevator position jumps between fixed values 
approximately every 10 sec. These jumps are the characteristic response expected for the 
modeled hysteresis where the elevator position only moves when the elevator command exceeds 
the amount of the hysteresis. To confirm that the stable low-amplitude oscillation was caused by 
the elevator acturator hysteresis, the simulation run was repeated with the elevator hysteresis set 
to zero. The results of this run are plotted in Figure 32 and these results show no oscillation in 
the altitude error. A method or technique to remove this oscillation was not found as part of this 
control law development effort. Note that the altitude error oscillation would not be seen by the 
pilots since the amplitude of the altitude error is below the resolution of the altitude displayed in 
the cockpit. It is possible that the stable oscillation may not be perceptible to pilots or passengers 
in the presence of normal environmental disturbances given the low magnitude and long time 
period of the oscillation. The maximum amplitude of the change in vertical acceleration during 
the oscillation was 0.003 g. 
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Figure 29. Various longitudinal variables for 37,000 ft cruise altitude capture. 
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Figure 30. Path errors and command signals for 37,000 ft cruise altitude capture. 
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Figure 31. Logic discretes for 37,000 ft cruise altitude capture. 
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Figure 32. Path errors and command signals for 37,000 ft cruise altitude capture with no elevator actuator 
hysteresis. 
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Speed Control Performance 


This report shows the results of four simulation runs to illustrate the performance of the speed 
control mode of the VNAV control system. The runs included are (1) CAS commanded to 311 
knots from 250 while initially climbing from 11,000 ft, (2) repeat of run 1 but ending with a 
constraint altitude capture, (3) a climb from level flight at 10,000 ft with CAS commanded to 
328 knots from 255 knots, and (4) an idle TOD run from 37,000 ft following a Mach/CAS 
schedule. The simulation data for these runs is shown as plotted data, and the same variables and 
discretes that were plotted to illustrate path mode control performance were plotted for 
illustration of the speed control mode performance. 

Figures 33-35 show the simulation data for the first run. The simulation begins with the aircraft 
trimmed at an altitude of 1 1,000 ft, CAS of 250 knots, flight path angle of 3.2124 deg, weight of 
195,000 lb, and a CG located at 19 percent MAC. At 0.22 sec, the VNAV control system is 
anned, current path segment 13 of the VNAV trajectory shown in Table 7 is captured, and the 
control system is engaged in the path control mode as shown in Figure 35 and begins tracking 
current path segment 13. The autothrottle is initially engaged in the speed control mode (Figure 
35) (note SPD_R). At 7.5 sec, the next segment path capture discrete is triggered and that 
causes the desired control mode and commands to be taken from path segment 14. Path segment 
14 requests the VNAV speed mode to be engaged, the CAS command to be 311 knots, and the 
autothrottle to be set to the constant thrust mode. As a result, the path mode disengages, and the 
VNAV speed mode is engaged. Also, the autothrottle is set to the thrust mode and that causes 
the engine PLA to be commanded to 106.15 deg for enabling the aircraft to both climb and 
accelerate (Figure 33). The vertical speed submode also engages at this time (Figure 35) since 
the true airspeed error exceeds the value specified by parameter DelTASVS (set to 42.195 ft/sec 
or 25 knots). When climbing, the vertical speed is commanded to the current vertical speed with 
a limit of 500 ft/min. For this run, since the vertical speed at speed mode engagement was 
greater than 500 ft/min, the vertical speed is commanded to 500 ft/min. The aircraft pitches 
down to control the vertical speed to that limited commanded value accelerating toward the 
commanded airspeed. The plot of vertical speed (Hdot) in Figure 33 shows that the vertical 
speed stays slightly above the dashed straight line in the plot which is equivalent to 500 ft/min. 
The vertical speed does not reach the 500 ft/min command because the aircraft is accelerating 
which requires the aircraft to constantly pitch down to maintain a constant vertical velocity. 
When the true airspeed error becomes less than 42.195 ft/s, the vertical speed submode 
disengages and the speed-on-elevator mode computes pitch steering commands to follow the 
commanded airspeed. Figure 35 shows that the vertical speed submode disengages around 38 
sec. At that time, the aircraft pitches up to reduce the flight path acceleration and capture the 
commanded CAS (Figure 33). Figure 34 shows that the capture of the commanded CAS is 
smooth with a 1 .2 knot overshoot around 65 seconds and by 80 sec the magnitude of CAS error 
is reduced to a value of about 0.5 knots. Note that for all simulation runs, the vertical 
acceleration parameter, SpdHddLim, was set to 3.0 ft/s 2 . 
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Figure 33. Various longitudinal variables for CAS command from 250 to 311 knots in speed control mode. 
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Figure 34. Path and speed errors and command signals for CAS command from 250 to 311 knots in speed 

control mode. 
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Figure 35. Discretes for CAS command from 250 to 311 knots for speed control mode of VNAV control 

system. 
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The second simulation illustrates the VNAV control system performance of a constraint altitude 
capture when the VNAV control system is engaged in the speed control mode. The run is a 
repeat of the run above except that a constrained altitude capture occurs around 53 sec. This run 
demonstrates constraint altitude capture which is a path control mode; however, the simulation 
run is included in this section because the altitude capture is engaged from the speed control 
mode and begins the same as the prior simulation run. Figures 36-38 show the results for this 
run. The discretes of Figure 38 show that the constrained altitude capture occurs at 53 sec when 
the aircraft is climbing with a vertical speed of nearly 34 ft/s (2,040 ft/min) and with a vertical 
acceleration of 2.5 ft/s 2 (Figure 36). The thrust control mode changes from fixed thrust to speed 
control at constraint altitude capture and the commanded PLA position begins to decrease to 
control speed to the 311-knot command. The control system begins to pitch the aircraft down 
when the altitude capture is engaged, but the vertical speed reaches nearly 38 ft/s (2,280 ft/min) 
before it begins to decrease. The vertical acceleration reaches a value of -3.0 ft/s 2 at 58 sec 
which is nearly twice the desired limit of 1.6 ft/s 2 (PthHddLim was set to 1.6 ft/s 2 for this run). 
The vertical acceleration limit is exceeded because, most likely, the path capture criterion was 
derived on the assumption that vertical acceleration would be near zero when path segment 
capture initiation occurred. At path segment capture initiation, the vertical acceleration was 2.5 
ft/s 2 . The plot of altitude in Figure 36 shows that the altitude capture is smooth, and the plot of 
altitude error in Figure 37 shows that the altitude command is overshot by 8 ft at 100 sec and 
probably would have reached about 10 ft before beginning to converge on the commanded value 
if the simultion time were longer. 
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Figure 36. Various iongitudai variables for constraint altitude capture during speed control mode 

engagement. 
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Figure 37. Path and speed errors and command signals for constraint altitude capture during speed control 

mode engagement. 
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Figure 38. Discretes for constraint altitude capture during speed control mode engagement. 
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The third simulation run illustrates the VNAV speed control mode performance in capturing a 
commanded CAS increase when starting from level flight at 10,000 ft with the throttle 
commanded to increase to a fixed position to begin a climb. For this run the aircraft is trimmed 
at an altitude of 10,000 ft, flight path angle of zero degree, CAS at 254.82 knots, weight 198,968 
lb, and a CG located at 20 percent MAC. The aircraft initial position is located at the beginning 
of path segment 10 for the VNAV trajectory in Table 6. Figures 39-40 show the results for this 
run for which the VNAV control system is armed for engagement at 0.22 sec, and at that time 
engages in the speed control mode and the vertical speed submode (Figure 41). Also, at that 
time, the throttle is commanded to a fixed throttle position of 104.09 deg and CAS is 
commanded to 328 knots (Figure 39). Figure 39 shows that the CAS begins to increase as a 
result of the throttle increase and the vertical velocity initially increases because of the pitchup 
caused by the thrust increase. Then, the vertical velocity reduces toward the commanded value 
of zero. The commanded value is zero because the aircraft vertical speed is zero when the 

vertical speed submode engages. At 31 sec, the vertical speed submode disengages when the 

true airspeed error becomes less than 25 knots and from that time the speed control law 
commands pitch steering signals to capture and track the commanded CAS. Figure 39 shows 
that the capture of the commanded airspeed is smooth, and the plot of CAS error in Figure 40 
shows that there is a 1.5 knot overshoot of the commanded CAS at 56 sec followed by a slow 
convergence towards zero error. At about 86 sec, the CAS error reaches 0 ft/s and then 
subsequently reaches a steady value of about 0.25 knots less than the commanded CAS. This 
small offset in this run is caused by the fact that, when the aircraft is climbing at constant CAS, 
true airspeed must be constantly increasing to maintain constant CAS. Note that the plot of 
flight path acceleration in Figure 39 shows that after constant CAS is established, it has a small 
positive value (a little less than 0.5 ft/s 2 ) which means that true airspeed is increasing. Thus, for 

a fixed thrust, the aircraft must pitch down to increase true airspeed. Also, in climbing at a fixed 

throttle position, the engine thrust is decreasing and a further reduction in pitch angle is needed 
to maintain true airspeed. So, the small lag in the aircraft pitch down causes the CAS to be offset 
from the commanded value. Note that without the use of the bias term in the speed command 
processor, described earlier in the Speed Control Law section of this report, the CAS offset 
would be larger. 
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Figure 39. Various iongitudai variables for climb from 10,000 ft and CAS increase in speed control mode. 
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Figure 40. Path and speed errors and command signals for climb from 10,000 ft and CAS increase in speed 

control mode. 
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Figure 41. Discretes for climb from 10,000 ft and CAS increase in speed control mode. 
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The simulation results for the fourth run in Figures 43-44 illustrate the performance of control 
system speed control mode for an idle descent from a cruise altitude. For this run, the aircraft 
model is trimmed at an altitude of 37,000 ft in level flight with a Mach number 0.7964, weight of 
172,841 lb, and a CG located at 25 percent MAC. The initial aircraft position is at the start of 
path segment 1 for the VNAV trajectory in Table 5. The VNAV control system is armed for 
engagement at 0.22 sec and, at that time, immediately engages in the speed control mode to track 
the commanded Mach number of 0.7964. Also, at this time, the throttle is commanded to idle at 
-2 deg/s. The plot of the commanded Mach minus actual Mach in Figure 43 shows that the 
Mach number decreases by 0.007 from the commanded Mach during the pitch down from 0 to 
about 25 sec with the throttle decreasing to idle. This error could be reduced by decreasing the 
magnitude of the throttle reduction rate or pitching down at a higher rate. Reducing the throttle 
rate is probably the best solution, since the pitch down rate resulted in a 0.1 g maneuver, as 
shown by the plot of vertical acceleration in Figure 42. For this run, the vertical acceleration 
parameter SpdHddLim for the speed control law is set to 3.0 ft/sec 2 . After the initial pitch down, 
the speed control system reduces the Mach error to less than 0.001. At 146 sec, MACHSEL 
becomes false, meaning that now the CAS command of 300 knots is to be tracked. The plot of 
CAS in Figure 42 and the plot of CAS error in Figure 43 show that there is a 2-knot overshoot in 
capturing the commanded CAS. This overshoot occurs because the change from tracking a 
constant Mach number to tracking a constant CAS requires a change in the flight path 
acceleration of the aircraft. The plot of flight path acceleration in Figure 42 shows that for a 
descending flight path at constant Mach number, the aircraft has a small acceleration and, when 
tracking constant CAS, has a small deceleration (about twice the magnitude of that for constant 
Mach number). So, if the MACHSEL discrete were set to false a little earlier by some 
anticipation calculations and logic, the CAS overshoot could be reduced. For zero CAS 
overshoot the flight path acceleration would have to change instantaneously from the positive 
value for constant Mach tracking to the negative value for constant CAS tracking. Currently, for 
this simulation, the MACHSEL discrete changed at the point where the true airspeed value for 
the commanded Mach and the true airspeed value for the commanded CAS are the same. The 
plot of CAS error in Figure 43 shows that after the initial CAS overshoot, the speed control mode 
maintains the magnitude of the CAS error at less than 0.25 knots. 
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Figure 42. Various longitudinal variables for idle TOD in speed control mode. 
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Figure 43. Path and speed errors and command signals for idle TOD in speed control mode. 
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Figure 44. Discretes for idle TOD for speed control mode. 
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Summary 


An outer loop VNAV control system was designed and developed to track VNAV trajectories 
intended for use in high-density traffic operations of the NextGen air traffic management system. 
This VNAV control system has path and speed modes. The speed mode has a vertical speed 
submode that engages when the speed error is large and counter in energy to the specified flight 
phase (i.e. accelerating while climbing or decelerating while descending). The VNAV control 
system also contains control logic to arm and engage the control modes and to initiate capture of 
VNAV trajectory path segments. The VNAV control system also contains control logic to level 
an aircraft using the path control mode to the constraint altitude set on the MCP if the aircraft is 
about to fly through the MCP altitude while tracking a VNAV trajectory. Therefore, the VNAV 
control system handles all nonnal VNAV operating modes. 

The perfonnance of the VNAV control system was shown using a Simulink six DOF, non-linear 
aircraft simulation containing a proprietary model of a B757-200 aircraft. The performance of 
the path control mode was shown for four test scenarios, and the performance of the speed 
control mode was shown for four different test scenarios. The test scenarios used trajectories 
generated for NextGen operations. Path tracking altitude errors of 5 ft or less were shown after 
path capture maneuvers and transients because of flap deployment. Speed tracking errors of 0.5 
knots or less were shown after capture of speed command changes. The simulation scenarios 
also showed transient-free engagement of the control modes for VNAV control system 
engagement, transition between the path and speed control modes, and transition to and from the 
vertical speed submode of the speed control mode. 

A stable low-amplitude oscillation on altitude error with an amplitude of ±1 ft was encountered 
on one test scenario after cruise altitude capture. The oscillation was determine to be caused by 
hysteresis modeled in the elevator actuator. A solution was not determined as part of this control 
law development effort. Note that the magnitude of the altitude error would not be be observed 
by the pilots since the resolution of the altitude displayed in the cockpit is an order of magnitude 
larger than the altitude error of this oscillation. Further, this oscillation may not be sensed by the 
pilots or passengers under normal environmental disturbances since the maximum change in the 
magnitude of vertical acceleration during the oscillation was 0.003 g. 

A stability analysis of the VNAV control system was conducted to determine if the gain and 
phase margins for the path and speed control modes satisfied nominal gain and phase margin 
requirements. For the path control mode, the minimum gain margin was 12 dB which exceeded 
the gain margin requirement of 6 dB and the minimum phase margin was 80 deg that also 
exceeded the phase margin requirement of 30 deg. For the speed control mode, the minimum 
gain margin for the sampled operating conditions was 10 dB and the minimum phase margin was 
80 deg. Since the gain and phase margins exceeded the gain and phase margin requirements by a 
fair amount, the VNAV gains could be increased to potentially improve path and speed tracking 
performance if so desired. However, gain increases would result in more aggressive path and 
speed captures with vertical accelerations expected to be larger. Since vertical acceleration 
exceeded desired limits in some of the performance results with the current gains, some 
additional means for maintaining vertical acceleration within desired limits would probably be 
required for increased gains. 
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Appendix A: Derivation of Estimated Delta Thrust to Weight Parameter 


The longitudinal equation of motion for a point-mass aircraft model is given as 

mV = T -D-Wsiny (Al) 


where m is the mass of the aircraft, V is the time derivative of the velocity vector, T is the 
aircraft thrust, D the aircraft drag, W the aircraft weight, and y is the flight path angle. The 
vertical velocity of the aircraft is given as 

h = V sin y (A2) 


where V is the total velocity vector. The mass of the aircraft, m, can be expressed as W/g where 
g is the acceleration due to gravity. Solving for sinyin eq. A2, substituting for that result in eq. 
Al, substituting W/g for in in eq. Al, and rearranging give the equation for the estimated delta 
thrust to weight parameter EstDTW as 


EstDTW = 


(T-D) 

w 


V h 
— + — 


8 V 


(A3) 


Note that multiplying EstDTW by g gives 


EstDTW g= V + h^ (A4) 

where the first term on the right of eq. A4 is the aircraft’s flight path acceleration, and the second 
tenn can be considered an estimate of an additional amount of potential flight path acceleration 
of the aircraft. The second term is considered potential acceleration since pitching to reduce the 
magnitude of the vertical velocity would result in a positive flight path acceleration for positive 
vertical velocity and negative flight path acceleration for negative vertical velocity. 


69 



Appendix B: Derivation of Rate of Change of True Airspeed for Constant 
CAS and Constant Mach Conditions 

The time rate of change of true airspeed as a function of altitude rate was desired for use in the 
second order speed command processor for the speed-on-elevator control law. The rate of 
change of true airspeed was desired for aircraft operations under either constant calibrated 
airspeed or constant Mach operating conditions. The input to the second order processor is the 
desired true airspeed that is detennined from either commanded calibrated airspeed (CAS) or 
commanded Mach number. The commanded CAS or Mach number is set either by a pilot via 
the aircraft MCP or by the flight management system. When the aircraft climbs or descends 
under a Mach hold or CAS hold condition, the true airspeed changes for all constant CAS or 
Mach conditions except for constant Mach above 36,089 ft. The Mach number above 36,089 ft 
remains constant since the atmospheric temperature is constant per the 1976 standard 
atmospheric model (ref. 7). 

When true airspeed is changing during a CAS or Mach hold condition, the output of the second 
order processor will lag the commanded input without compensation. Computing the rate of 
change of true airspeed as a function of altitude rate and inserting the computed value at the 
input of the output integrator of the second order processor can overcome this lag. The 
computation of true airspeed as function of altitude rate follows. 


True airspeed is defined as 


V T = aM = a 


7-1 


r - 1 


^ + 1 


-1 


(Bl) 


where a is the speed of sound at some atmospheric model altitude, M is the Mach number, 7 is 
the ratio of specific heats (c>/c\), q, is impact pressure, and P is static pressure at the atmospheric 
model altitude. For air, 7 is equal to 1.4. 

Derivations for Constant Calibrated Airspeed 


Calibrated airspeed is defined as 


V=a 


std i 


7-1 


r- 1 


-^ + 1 

P 

std J 


-1 


(B2) 


where astd is the speed of sound at sea level standard conditions and Pstd is the static pressure at 
sea level standard conditions. Note that the speed of sound is given as 


a = a 


std 



(B3) 


where T is the temperature at atmospheric model altitude and Tstd is the atmospheric temperature 
at sea level standard conditions (Tstd =518° Rankin). 
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From eq. Bl, the change of true airspeed with altitude, dVr/dh, is 


Let 


and 


dV, dM da 

— - = a + M — 

dh dh dh 



c 2 


7-1 

r 


so that 


M = 



(B4) 


(B5) 


(B6) 


(B7) 


Change in Mach Number with Altitude for Altitudes Equal to or Below 36,089 ft 


The change in Mach number with altitude, dM/dh, under the conditions for constant calibrated 
airspeed for altitudes of 36,089 ft or less can be found via eq. B7. Note that for constant 
calibrated airspeed, q c must be constant (eq. B2). Differentiating eq.B7 with respect to altitude 
gives 


-i 


dM 1 

f ■ 

f C P ^ 

c 2 

\ 

2 

f q c 

= — 

c. 



-1 

CaCj 


dh 2 

l 

l . 

[P J 


J 

1 z 

l P J 


(-Ik ?- 2 


dP_ 

dh 


(B8) 


or with simplification yields 


dM 

dh 


C 1 C 2 

2 M 


f 

v 



Y ' 2 -1 

J 



dP 

dh 


(B9) 


Now, derive dP/dh. For the 1976 atmospheric model we have the relationship between the static 
pressure ratio and temperature ratio for altitudes equal to or less than 36,089 ft as 



/ 

V 


b?R 


1 std J 


(BIO) 


where 


T = bi + bah 


hi = 518° Rankin 
bi = -0.003566 deg Rankin/ft 
R = 53.35 ft/deg Rankin 
h = altitude, ft 
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Thus, upon substitution 


P = P„ 


( r r ^ 


\ T stdJ 


std 


( b, + b 2 h ^ 

V T std J 


where 

b3=-\lb 2 R = 5.256 
Given eq. B1 1, dP/dh can be found as 


dP h A P s t d 

( b, + b 2 h } 

i 3 -1 

b h P 

_ u 2 U Z r std 

( j ^ 

dh T std 

v T std J 

T 

std 

l P sul y 




Solving for T from eq. B1 1 gives 


T-T 


std 


f _P_ A 

V J 


Substituting for T in eq. B12 gives 

f 1 V’a -1 

~i 

P 

= - J -“ | | ■ 

dh 


dP b b P 

ur _ v 2 u 3 r std 


' S/rf 


Shi 


V 


h h P 

u 2 u 3 r std 


■ sfc/ 


y 


A D A 


V P s,d J 


b 3 - 1 


1/i, 




h h P 

u 2 u Z r std p b 3 


std 


Now slightly rearrange eq. B9 to give 


dM -cx 


W2 


dh 2M 


^ + 1 


Kp-idP 

P dh 


Now solve for qdP in terms of M. Square both sides of eq. B7 to give 


1 

f ■ 

+ 

o| 

<=2 ^ 

II 

CM 

c. 

-1 

1 

1 

l . 

\P ) 

-) 


After manipulation we get 


^ + 1 = 


^ M 2 


V c i 


A 


+ 1 


Substituting for qdP + 1 and qdP in eq. B15 gives 
dM 


" C 1 C 2 



c 2 -1 



+ 

CM 

c 2 


S' 

+ 

CM 

'v C 1 ) 



l c i ) 


dh 2 M 
Substituting for dP/dh from eq. B14 gives 


p - 1 dP 

dh 


dM —c^c 2 b 2 b 2 P s 


1/*3 

std 


dh 


2MT 


std 



c 2 



f M2 

c 2 


( M 2 P 

+ 1 



+ 1 

l c i J 



v c i y 


-1 


p-1lb 3 


(Bll) 

(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 

(B18) 

(B19) 
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From eq. B13 we get 




std 


( T \>*T Vh 

V ^ std J 


p-1/^3 
^ std 


( T V 1 

V T std j 


Substituting the speed of sound, a, from eq. B3 into eq. B1 gives 


^ T~ a stdJ T M 

V 1 std 


Squaring both sides of eq. B21 above and solving for T gives 


T = 


T V 

1 std * T 

al d M 2 


Substituting for T in eq. B20 gives 


p-1/^3 _ p-'Vb 3 
1 r std 


r y2 V 1 
K a st d M 2 j 


p -MH a ld M2 

*std T t 2 


v; 


P -MH a lM 2 
std V T V T 


Substituting the definition of Vt from eq. B1 for one Vt into the equation above gives 


p~Mb 3 _ p-Vb 3 a ld M 

1 r std x r 

V T a 


Substituting for P vt>3 in eq. B19 gives 


dM -c^c 2 b 2 b 3 a 


dh 


std 


2T s,d aV T 



C 2 -1 

~ 


f M 2 

c 2 


f M 2 h! 

+ 1 



+ 1 

V c i J 



V C 1 J 


After multiplying terms in the bracket by first term in parentheses of eq. B25, we get 


dM -c 1 c 2 h 2 h 3 a 


dh 


2T„ d aV r 


c 2 - 1 1 


[ M 2 } 

c 2 c 2 

{ M2 i) 

+ 1 

~ 

— + 1 

V C 1 J 


V C 1 J 


c 2 "1 

^2 


or 


dM _ -cY 2 b 2 b^a std 


dh 2T sld aV r 


+ 

CM 


f M 2 

+ 1 

V c i ) 


l c i J 


c 2~1 

c 2 


(B20) 

(B21) 

(B22) 

(B23) 

(B24) 

(B25) 


(B26) 
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The second term in parentheses may be converted to a series by the binomial theorem that is 
expressed as 


i (n - 1 )(n - 2)x 3 


2! ^ 3! 

n{n - \)(n - 2) - ■ -{n - r + \) x r 

rl 


This series converges if \x\ < 1, and since \M 2 /ci\ is less than one, the series for the second tenn in 
parentheses of eq. B28 will be a convergent series and is given as 




(B28) 


(B29) 


The static pressure, P, for altitudes above 36,089 ft is given by 


P T 

p _ r std 1 36 


- 0 - 36 , 089 ) 


where Tu is the atmospheric temperature at 36,089 ft (390° Rankin), R is the gas constant for air 
(53.35 ft/° Rankin), h is the atmospheric altitude, and e is the base of the natural logarithm. 
Differentiating eq. B30 with respect to h gives the pressure change with altitude as 
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m -(*-36,089) 

dP _ Pstdl 36 RT m 

{ _1 ] 

= P 

r -i ^ 

dh T std 

V ^ ^ 36 J 

V ^^36 y 


(B31) 


Now substituting dP/dh above into eq. B18 gives 


dM -qc 2 


dh 2 M 


M 


V c i 


+ 1 


r M 2 


V c i 


+ 1 


-1 


f -1 A 


V ^^36 J 


(B32) 


After some manipulation of eq. B32, use of the binomial theorem, and neglecting higher order 
tenns as was done to get eq. B29, eq. B32 can be expressed as 


dM 


1 


dh 2RT x M 


m 2 + 


r c 2 ~^ 


2qc : 


M 4 + 


2 J 


-cp 

6 c 2 c\ j 


M 


(B33) 


Since the atmospheric temperature remains constant above 36,089 ft, the speed of sound remains 
constant and can be defined as 


^36 — °std i 


' 36 


1 Std 


(B34) 


and da/dh is zero above that altitude. Substituting for M equal to Vrlcm in the denominator of eq. 
B33 gives 


dM 


36 


dh 2RT 3e V r 


M 2 + 


f c 2 -^ 


M 4 + 


2 qc 2 j 

Derivations for Constant Mach Number 

For constant Mach number, eq. B4 becomes 

dV T da 

— - = M — 
dh dh 


-cp 

6c 2 c 2 j 




(B35) 


(B36) 


Thus, a solution requires the derivation of da/dh. The speed of sound can be expressed as 

a = ^ygRT =jygR(b,+b 2 h) (B37) 

where /is the ratio of specific heats (cp/cv equal to 1.4), g is gravity constant ( g = 32.174), and R 
is the gas constant for air ( R = 53.35 ft/° Rankin), and T is the atmospheric temperature. 


Change in Speed of Sound with Altitude for Altitudes Equal to or Below 36,089 ft 
Differentiating eq. B37 with respect to h gives 


(B38) 
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(B39) 


Substituting for a as defined by eq. B1 gives 
da ygRb 2 M 
dh ~ 2 V t 

Change in Speed of Sound with Altitude for Altitudes Above 36,089 ft 

As noted earlier, above 36,089 ft temperature, and thus the speed of sound, remains constant and 
da/dh = 0. 

Time Rate of Change of True Airspeed for Constant Calibrated Airspeed and Constant 
Mach Number Conditions 


The time rate of change of true airspeed, dVr/dt, can be derived by substituting the derived 
expressions for dM/dh and da/dh into eq. B4. 

dV r /dt for Constant Calibrated Airspeed 


The time rate of change of true airspeed can be written as 
dV T dV T dh dV T ^ 

dt dh dt dh 

The equation for dVr/dh is repeated here from eq. B4 

dV T dM da 

— -- a + M — 

dh dh dh 


Substituting for dVr/dh in eq. B40 gives 


dV T 

dt 


dM 

> 

dh 


da 
■ + M — 


dh 


h 


(B40) 


(B41) 


Altitudes Equal to or Less Than 36,089 ft 


Substituting dM/dh and da/dh from eq. B29 and B39 gives 


dV T 

dt 


h 2 b 3 a std 


IT 


std 


M 2 + 


r c 2 -0 
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M 4 + 


2 J 


Tzff 
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or combining terms gives 


dV T 

dt 


~b 2 b 3 a- std + ygRb 2 


9 T 2 

lL1 std z ' 

—b 2 b 2 a std ( 1 ~ c 2 j 

^std C 1 C 2 


M 2 + 


-bib/Cd (A ~ l) 

4T s,d C l C 2 


M 




h_ 

K 


(B42) 


(B43) 
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Altitudes Greater Than 36,089 ft 


Above 36,089 ft, the speed of sound is constant and, thus, da/dh is zero. Therefore, dVr/dt is 
given as 


dV T 

dt 


dM • 

h 

dh 


dM • 

h 

dh 


(B44) 


Substituting for dM/dh from eq. B35 gives 


dV T 

dt 


*36 


2 RT 


M 2 + 


36 


2cjC 


M 4 + 


2 J 


'izfT 

6c\c\ J 


M 


Wr 


dVj/dT For Constant Mach Number 


(B45) 


There are two conditions for determining the rate of change for true airspeed when changing 
altitude at constant Mach number - altitudes equal to or less than 36,089 ft and altitudes greater 
than 36,089 ft. 


Altitudes Equal to or Less Than 36,089 ft 


For constant Mach number, dM/dh is zero and, thus, from eq. B41 dVr/dt is given as 


dV T 

dt 


da 
M — 
dh 


h 


Substituting for da/dh from eq. B39 in eq. B46 gives 


dV T ygRh 2 M 2 


dt 


< it} 

\V T J 


(B46) 


(B47) 


Altitudes Greater Than 36,089 

Since da/dh is zero for altitudes greater than 36,089 ft for constant Mach, dVr/dt is also zero for 
these altitudes. 


M 2 , M 4 , and M 6 Coefficients 

Let the definition of dVr/dt for constant CAS at altitudes equal to or less than 36,089 ft be 

^ = (k l M 2 +k 2 M 4 +k 6 M 6 )— (B48) 

dt Vj 

where ki, ki, and h are constant coefficients. 

Let the definition of dVr/dt for constant CAS at altitudes greater than 36,089 ft be 


77 



(B49) 


^ = (k,M 2 +k,M 4 +LM 6 ) — 

clt V 3 4 7 ’v T 

where ks, k4, and k? are constant coefficients. 

Let the definition of dVr/dt for constant Mach at altitudes equal to or less than 36,089 ft be 

^Yl = U 5 M 2 )— (B50) 

dt v 5 ’v T 

where ks is a constant coefficient. 

Then, the coefficients for the M 2 , M 4 , and M 6 terms in eq B48-50 above can be determined from 
eq. B43, B45, and B47. These coefficients are summarized in the table below. 


Table 8. Constant coefficients for computing dVj/dt for constant CAS and Mach conditions. 


Terms 

Constant CAS 

Constant Mach Number 


h < 36,089 ft 

h > 36,089 ft 

h < 36,089 ft 

h > 36,089 ft 

M 2 

, _ ~ b 2 b 3 a ii , rs Rb 2 
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r _ ^36 
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2RT X 
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4 47CT c c 

36 c 1 c 2 
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0 

“ M * 

~ b 2 b 3 a ld 0 - C 2 ) 
6 qt 2 2 

®T std c 1 C 2 
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, ~ C 2 ) 
7 ' 8 RT m c\c\ 

k? = 2.5371 

0 
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Note that substituting for bs as defined in this appendix after eq. BIO and au from eq. B34, the 
constant coefficients for constant CAS become 


! 7***2 
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Thus, it is observed that ki = k.4 and k6 = h 



Appendix C: Autothrottle Used for VNAV Control Law Performance Testing 

Figure 45 shows the autothrottle control law used for VNAV control law performance testing. 
This control law is a modified version of the control law that was developed for the research 
flight deck of the B737 aircraft described in ref. 1. The output of the control law is the 
commanded position of the power lever angle (PLACMD) for input to the B757 autothrottle 
servo model in the Simulink simulation. The control law either tracks the commanded Mach 
number (MachCmd) or commanded calibrated airspeed (CASCMD) when the autothrottle is 
controlling the airspeed (SPDR is true). When, SPDR is false, the commanded PLA will be 
constant for emulating a fixed EPR condition or be driven to a commanded idle position 
(PLAIDLE) at a rate of -2 deg/sec when IDLE R becomes true. 



i.c. = PLA 

@ IDLE_R • VNAVE 


Figure 45. Autothrottle of B737 research airplane. 
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Appendix D: VNAV Simulink Block Diagrams 


The VNAV control system was developed using an existing 6 -DOF non-linear aircraft Simulink 
simulation where the aircraft dynamics represent the Boeing 757-200 aircraft with Rolls Royce 
engines. The VNAV control system was developed in the Simulink environment of this 
simulation. The simulation includes non-linear models of B757-200 autopilot servos, surface 
electro-hydraulic actuators, autothrottle servos, and the Rolls Royce engines. The description of 
Simulink diagrams and Simulink VNAV control system that follows assumes the reader is 
familiar with the S im ulink software environment. Figure 46 shows the top level block of the 
B757 Simulink simulation. The VNAV control system is implemented in the block labeled 
Autopilot or FCC at the lower left of the top level block. This block is provided for autopilot 
control law implementation in this simulation. 


©■ 





Cockpit Switches 
and Misc. Values 



Autopilot 

Engage 



Simulation 

Options 


Autopilot or FCC 


Figure 46. Top level Simulink block of B757 Simulink simulation. 


Figure 47 shows the contents of the block Autopilot or FCC. The VNAV control system is 
contained within the blocks PathGuid, VNAV_Mode_Logic, VNAV_Control_Law, and 
GainScheduleParams. Some of the remaining Simulink blocks contain the B757 pitch and roll 
inner loop control laws and computed parameters for these control laws. The details of these 
blocks are proprietary to Boeing and are not described here. The block labeled FlapCommand 
was implemented to command flap extension or retraction based on CAS landing reference 
speeds for the B757. The block labeled AutoStab contains the B757 autotrim algorithm that 
commands the trim actuator motors to move the stabilizer surface for either a trim up or trim 
down deflection to keep steady-state deflection of the elevator surface close to the neutral 
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position. The block labeled PitchCmdProcess (a Boeing proprietary block) computes a pitch 
command from the VNAV control system pitch steering command signals VNAVS and VNAVI. 
The block labled PulseLinDisc is used to interface the autopilot surface commands to the 
autopilot servo actuator models and the auto stabilizer trim commands to the model of the 
stabilizer trim motors. The autothrottle used for the VNAV control law development is 
contained in the block labeled AutoThrottle. 



Figure 47. Contents of Autopilot or FCC Simulink block. 


Figure 48 shows the contents of the block PathGuid. The blocks labeled PathCaptureLogic and 
ErrorSwitch are part of the VNAV control system. The blocks could be considered part of the 
RPFMS calculations, but would require passing different variables and discretes between the 
RPFMS calculations and the VNAV control system. The remaining blocks were developed to 
emulate the calculation of VNAV path errors and discretes that would be computed by the 
RPFMS VNAV guidance so that the VNAV control system performance could be determined 
under realistic scenarios. 
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— ► rwCoords 


Commands ~ 


— ► Sensors 


PathErrors 
EnableDiscs 
EnggeDiscs 
ModeReqsl 
PathCaptureLogic 


G> 

EngageDiscs 


CntlErrors 

PathStalus 

EnableDiscs 

EngageDiscs 

SwErrDiscs 

VNAV_EnableDiscs 


— ► SwErrDiscs 


— ► ModeParams 


— ► Commands 


VNAV_ModeDiscs 


EnableDiscs 

PathErrors CntlErrors 
PthCapDisc 
PathStalus 
ModeReqst 
SpdSel 

ErrorSwitch 


SwErrDiscs 


-*C D 

CntlErrors 


VNAVCntIDiscs 


>QD 

VNAVCntIDiscs 


— ► Sensors 


-K7) 

CMOS 


SpeedCalcs 


Figure 48. Contents of Simulink block PathGuid. 


Figure 49 shows the contents of block PathErrors in Figure 48. Block PathErrors contains a 
Matlab function block labeled Vertical Path Seg Calcs, which calls the Matlab function 
patherrors, and the block TrajectSigDef. Figure 49 lists the Matlab function patherrors. 

The Matlab function patherrors (1) computes the altitude and altitude rate errors for the current 
VNAV path segment, the next path segment, and the constraint altitude, (2) sets the desired 
control mode parameters and speed commands for the path segments, and (3) determines the 
path status. The details of these four items are presented in the discussion of items that are 
contained in block TrajectSigDef. The inputs to the function are the VNAV trajectory data 
(pathdata) that contains data such as that shown in Tables 5, 6, and 7, the past value of the path 
segment number (isegp), the ground range relative to the path segments (xa), aircraft altitude 
(alt), aircraft vertical speed (hdot), aircraft ground speed (vg), and the constraint altitude (altsel). 
The Simulink constant block MCPaltsel provides for setting the constraint altitude for a 
simulation run. The constant block Range provides a means for setting the aircraft ground range 
relative to the VNAV trajectory range data in the pathdata. For example, pathdata is a two- 
dimensional array where the first array index is the VNAV path segment number and the second 
index represents various information about that path segment (see Table 5). The third column 
contains the accumulated ground range at the start of the path segment. The memory block 
labeled PastSegNum is used to set the initial VNAV path segment number and to provide the 
past path segment number computed for input to the function patherrors. The outputs of the 
function patherrors will be addressed in the discussion of the contents of block TrajectSigDef. 
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MCPaltsel 


MCPaltsel must be set to the altitude flying to 
even if the control model is not the path mode 
because the computed limits in the 2nd-order 
processor for the speed mode change between 
flying up to an altitude or flying down to an altitude 


Figure 49. Contents of Simulink block patherrors. 

The Matlab function patherrors that is called in Simulink block Vertical Path Seg Calcs is listed 
below. 

function [ PathErrors ,ModeParams , Commands , PathStatus ] = 
patherrors (path, isegp, xa , alt , hdot , vg , altsel ) 

% 

% This block supports the Embedded MATLAB subset. 

% See the help menu for details. 


% 

% The variable path(i,j) contains information about the vertical path 
% segments where i represents the segment number and j is defined as 
% follows: 


% 

% 

% 

% 


% 

% 


path(iseg,l) = the Excel data file line number 
path(iseg,2) = tangent of the flight path angle of the ith 
path(iseg,3) = is the ground range of the start of the ith 
path(iseg,4) = altitude at the start of the ith segment 
path(iseg,5) = flight phase of the ith segment 
path(iseg,6) = pitch mode of the ith segment 
path(iseg,7) = thrust or throttle mode of the ith segment 
path(iseg,8) = CAS command of the ith segment 
path(iseg,9) = Mach command of the ith segment 
path ( iseg, 10 ) = flap command of the ith segment 
path ( iseg, 1 1 ) = gear command of the ith segment 


segment 

segment from the origin 


% The flight phase, pitch mode, and throttle mode are defined as follows: 
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% 

% 

% 

% 

% 

% 

% 

% 

% 

% 

% 


path ( iseg, 5 ) 
path ( iseg, 5 ) 
path ( iseg, 5 ) 


1 is CLIMB phase 

2 is CRUISE phase 

3 is DESCENT phase 


path ( iseg, 6 ) 
path ( iseg, 6 ) 


0 is VNAV speed mode 

1 is VNAV path mode 


path ( iseg, 7 ) 
path ( iseg, 7 ) 
path ( iseg, 7 ) 


1 is fixed thrust mode (e.g. climb thrust) 

2 is speed mode 

3 is IDLE mode 


% 

% 

% 

% 

% 

% 

% 

% 


xa - ground range from initial path waypoint, ft 

alt - altitude of aircraft, ft 

hdot - vertical speed of aircraft, ft/sec 

vg - ground speed of aircraft, ft/sec 

isegp - past value of path segment number 

nsegcapp - past value of next segment capture discrete 

asegcapp - past value of constrained altitude capture discrete 


J w LA kj L, V LA LA V LA La LA 11U V LA _L. 11L, LA LA V 

NAV Pth Eng p - past value of path mode 


% 

% If capture of the next segment occurred the past iteration and the 
% VNAV path mode is engaged, then update the segment number 

% 


PathErrors = zeros(6,l); 

ModeParams = zeros(6,l); 

Commands = zeros(8,l); 

PathStatus = zeros(3,l); 

iseg = isegp; 
segupd = 0; 

% Increment path segment is ground range is greater than ground range at the 
% beginning of the next path segment and not on last path segment 
if xa > path ( iseg+1 , 3 ) && isegp < ( size ( path , 1 ) -1) 
iseg = isegp + 1; 
segupd = 1 ; 

end 

%%%% set endseg to 1 if iseg value is last path segment in variable path 

if iseg == ( size (path, 1 ) - 1) 
endseg = 1 ; 

else 

endseg = 0; 

end 


% compute altitude error (DeltaAlt) and altitude rate error (DeltaHdot) 

% for current vertical path segment 

% 

he = path(iseg,4) + (xa - path ( iseg, 3 ) ) *path ( iseg, 2 ) ; % command altitude 

DelAltc = he - alt; 

DelHdotc = vg*path ( iseg, 2 ) - hdot; 

% 

% compute altitude error (DelAltn) and altitude rate error (DelHdotn) for 
% next vertical path segment 

% 

hen = path ( iseg+1 , 4 ) + (xa - path ( iseg+1 , 3 )) *path ( iseg+1 , 2 ) ; 

DelAltn = hen - alt; 

DelHdotn = vg*path ( iseg+1 , 2 ) - hdot; 
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% Compute altitude error (DelAltc) and altitue rate error (DelHdotc) for 
% constrained altitude mode 

DelAltcs = altsel - alt; 

DelHdotcs = -hdot; 

% Store current, next, and constrained altitude errors in output vector 

% 

PathErrors ( 1 ) = DelAltc; 

PathErrors ( 2 ) = DelHdotc; 

PathErrors ( 3 ) = DelAltn; 

PathErrors ( 4 ) = DelHdotn; 

PathErrors ( 5 ) = DelAltcs; 

PathErrors ( 6 ) = DelHdotcs; 

% 

% Set flight phase ( path ( iseg, 5 ) ) , pitch mode (path( iseg, 6 ) , and thrust 
% mode ( path ( iseg, 7 ) ) of the current & next segments 

% 


ModeParams ( 1 ) = path ( iseg, 5 ) ; 
ModeParams ( 2 ) = path ( iseg+1 , 5 ) ; 
ModeParams ( 3 ) = path ( iseg, 6 ) ; 
ModeParams ( 4 ) = path ( iseg+1 , 6 ) ; 
ModeParams ( 5 ) = path ( iseg, 7 ) ; 
ModeParams ( 6 ) = path ( iseg+1 , 7 ) ; 


% flight phase current segment 
% flight phase next segment 
% pitch mode current segment 
% pitch mode next segment 
% thrust mode of current segment 
% thrust mode of next segment 


% Store CAS command, Mach command. Flap command, & Gear command for current 
% and next path segments in output vector 


Commands ( 1 ) = path ( iseg, 8 ) ; 
Commands (2) = path ( iseg+1 , 8 ) ; 
Commands (3) = path ( iseg, 9 ) ; 
Commands (4) = path ( iseg+1 , 9 ) ; 
Commands (5) = path ( iseg, 10 ) ; 
Commands ( 6 ) = path ( iseg+1 , 10 ) ; 
Commands ( 7 ) = path ( iseg, 1 1 ) ; 
Commands ( 8 ) = path ( iseg+1 , 1 1 ) ; 


% CAS command of current segment 
% CAS command of next segment 
% Mach command of current segment 
% Mach command of next segment 
% Flap command of current segment 
% Flap command of next segment 
% Gear command of current segment 
% Gear command of next segment 


PathStatus ( 1 ) = segupd; 
PathStatus ( 2 ) = endseg; 
PathStatus ( 3 ) = iseg; 


Figure 50 shows the contents of TrajectSigDef. This block extracts the data from the output 
vectors of the function patherrors using Simulink demux blocks and puts that data into four buses 
using Simulink bus creator blocks. The figure shows that the bus PathError contains the altitude 
and altitude rate errors for the current and next path segments, and the constraint altitude along 
with the aircraft ground range. Aircraft ground range is included on this bus to provide later 
access for simulation data recording or to observe with Simulink tools such as Scope or Display 
blocks for debugging purposes. Simulink bus ModeParams contains the flight phase, desired 
pitch mode or VNAV control mode, and desired thrust or throttle mode for the current and next 
path segments. Simulink bus Commands contains the trajectory path segment commands for 
CAS, Mach number, flaps, and gear position for the current and next path segments along with 
the constraint altitude (AltSel). Note that these commands can be overridden by Simulink code 
in block VNAV ModeDiscs (to be discussed later in this section), and the flap commands can 
also be overridden by code within the Simulink block FlapCommand. The override capability 
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allows more flexibility in testing the VNAV control system. Simulink bus PathStatus contains 
the one-shot discrete segupd that signals when the aircraft reaches the starting point of a path 
segment, the discrete endseg that indicates if the aircraft is on the last path segment, and the 
current path segment number (iseg). 



Figure 50. Contents of Simulink block TrajectSigDef. 


Figure 5 1 shows the contents of block VNAV EnableDiscs. The block provides a means to 
emulate the setting of arming discretes, mode enabling discretes, and the VNAV disengage 
discrete. The discrete CurCapEnb enables capture of the current path segment when the 
VNAV control system is armed for engagement or armed for path mode engagement when the 
VNAV control system is engaged in the speed mode. The discrete NSEG_ARM enables the 
VNAV control system to capture the next path segment. The discrete Alt Cap Enb enables the 
VNAV control system to capture the constraint altitude, and the discrete EnableVSsubmode 
enables the VNAV speed control mode to engage the vertical speed submode. A value of one for 
these discretes represents the true state and value of zero the false state. 


87 



Figure 51. Contents of Simulink block VNAV EnableDiscs. 


Figure 52 shows the contents of block VNAV_ModeDiscs. This block allows for override of the 
control mode requests and commands from the VNAV trajectory data file. The block also sets 
the flight phase, requested control mode, and thrust mode when the constraint altitude is being 
captured or tracked (i.e. when CNTLASEG is true). Figure 53 shows the contents of block 
FLTPHASE which sets flight phase logical discretes from the flight phase integer variable 
defined by the VNAV trajectory data file. The flight phase integer from the VNAV trajectory 
file can be overriden with the switch and constant blocks shown at the top of the figure. Figure 
54 shows the contents of block PmodeRqst which sets logic discretes for requesting engagement 
of the path or speed control mode using the data defined by the integer variable in the VNAV 
trajectory data file. Figure 55 shows the contents of block ThrotMode which sets logic discretes 
for the desired thrust or autothrottle mode defined by the integer variable in the VNAV trajectory 
data file. 
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Pitchmode - VNAV control mode for current path segment - Path or Speed mode (0 - speed mode, 1 - path mode 
PitchmodeN - VNAV control mode for next path segment 
Fltphase - Climb, Cruise, Descent, or Altitude Hold flight phase 
VNAV_Diseng - Disengage VNAV Mode 

NSEG_ARM - A valid computation of the next segment path errors is being made 
VNAVARM - Arm VNAV mode for engagement 

MACHSEL - When 1 or true track the Mach command; when 0 track the CAS command 

SPD_R - Reguest autothrottle command throttles to track Mach or CAS 

THR_R - Reguest throttles be set at fixed thrust; e.g. max climb thrust 

IDLE_R - Request throttles be commanded to idle setting 



Figure 52. Contents of Simulink block VNAVModeDiscs. 



Figure 53. Contents of Simulink block FLTPHASE. 
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Figure 54. Contents of Simulink block PmodeRqst. 



Figure 55. Contents of Simulink block ThrotMode. 


Figure 56 shows the contents of block SpeedCalcs. This block emulates the RPFMS calculation 
of the selected true airspeed from selected CAS or selected Mach number. It also computes the 
MACHSEL discrete that determines if the selected CAS or selected Mach number is to be 
tracked by the control system. The block computes the true airspeed error DeltaTAS that is used 
by the VNAV speed control law to determine when to engage or disengage the vertical speed 
control submode. The From block, with variable name TASF, is the fdter true airspeed defined 
in Figure 4. The From block with variable name VTM is the computed conversion factor to 
calculate the true airspeed for the current Mach number measurement. The From block with 
variable name VTVC is the computed conversion factor to calculate the true airspeed for the 
current calibrated airspeed measurement. VTM and VTVC are computed in block 
GainScheduleParams in Figure 47. Block SpeedCalcs also computes the CAS and Mach errors 
(DeltaCAS & DeltaMach) that are used to record the performance of the VNAV speed control 
mode. 


90 



Figure 56. Contents of Simulink block SpeedCalcs. 

Figure 57 shows the contents of block PathCaptureLogic that is part of the VNAV control 
system. This block determines when a current path segment, next path segment, or constraint 
altitude is to be captured. How this logic works was discussed in the report section titled Path 
Capture Logic. The contents of the blocks One shot, One shotl, and One_shot2 are the same 
(Figure 58). The contents of the blocks KHERR1, KHERR2, and KHERR3 are all the same 
(Figure 59). 



Figure 57. Contents of Simulink block PathCaptureLogic. 
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Figure 58. Contents of Simulink block One shot. 



Figure 59. Contents of Simulink block KHERR. 


Figure 60 shows the contents of block Error Switch. This block is part of the VNAV control 
system and selects which set of altitude and altitude rate errors are used by the path control law. 
The errors that are selected are a function of the logic discretes CNTLNSEG and CNTLASEG. 
Figure 61 shows the contents of block CntlNsegCalc that compute CNTLNSEG.. Figure 62 
shows the contents of block CntlAsegCalc that compute CNTLASEG. A description of the logic 
is given in the report section titled Path Control Logic. 
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CntlAseqCalc 


Figure 60. Contents of Simulink block ErrorSwitch. 



Figure 61. Contents of Simulink block CntlNsegCalc. 
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Figure 62. Contents of Simulink block CntlAsegCalc. 


Figure 63 shows the contents of the block VNAV Mode Logic that is contained in block 
Autopilot or FCC (Figure 47). This block determines when the path control mode, speed control 
mode, or vertical speed control submode should be engaged or disengaged. The block also 
determines when to set the discrete PSYNC used for pitch steering command signal smoothing, 
integrator initialization, and rate limiter initialization. The calculations of this block are 
explained in the report subsection titled Engagement Logic. Figure 64 shows the contents of 
block EngageRequest, fig. 65 shows the contents of block VNAV_VS_Eng_Logic, and fig. 66 
shows the contents of block PSYNC Calc. 
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Figure 63. Contents of Simulink block VNAVModeLogic. 



Figure 64. Contents of Simulink block EngageRequest. 
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Figure 65. Contents of Simulink block VNAVVSEngageLogic. 



Figure 66. Contents of Simulink block PSYNC Calc. 
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Figure 67 shows the contents of the block Timer in Figure 66. VNAVD is set true 0.2 sec after 
VNAV engages and is used so that the signals labeled two, three, four, five and six in fig. 66 do 
not become true when VNAVE is first set true. 



Figure 67. Contents of Simulink block Timer. 


Figure 68 shows the contents of block VNAV Control Law that are contained in block 
Autopilot of FCC (Figure 47). The main parts of this block are the VNAV path control law in 
block VNAV_Pth ( Figure 69), the VNAV speed control law in block VNAV_Spd (Figure 71), 
and the pitch steering signals smoothing function in block FadelnOut (Figure 79). The 
functionality of these blocks is discussed in report section VNAV Control Laws. 



VNAV_Spd 


Figure 68. Contents of Simulink block VNAV Control Law. 
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Figure 69. Contents of Simulink block VNAV Pth. 


The contents of block KHERR2 in Figure 69 are the same as those for the KHERR1, KHERR2 
and KHERR3 for the path segment capture logic (Figure 59). Figure 70 shows the contents of 
the block Ratelim in the path control law block. The block Ratelim limits the input signal 
“dhlim” between rates ±PthHddFim, and the initial condition value (IC) is stored whenever the 
initial condition discrete (ICtrig) is set true. 



Figure 70. Contents of Simulink block Ratelim. 
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The contents of the blocks contained within block VNAV Spd are shown as follows: 
2ndOrderProcessor in Figure 72, VTerrorWashout in Figure 77, VertSpdControl in Figure 78, 
and Ratelim in Figure 70 (same contents of that Ratelim block in path control law of Figure 69). 



Figure 72. Contents of Simulink block 2ndOrderProcessor. 


The contents of the various blocks within the Simulink block 2ndOrderProcessor (Figure 72) are 
shown below as follows: block Limits in Figure 73, block Mid Value in Figure 74, block Bias 
in Figure 75 and block VarLimlnteg in Figure 76. As described earlier in the report sections, 
Speed Control and Second Order Speed Command Processor, the second order processor block 
smooths and limits the input command. The limits are acceleration limits whose values are those 
that the aircraft can physically attain. This smoothing and limiting is required since the true 
airspeed command input, TASsel, can have large step changes. The large changes in true 
airspeed commands can result from either changes in the pilot selection via the mode control 
panel or RPFMS command changes. 
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Figure 73. Contents of Simulink block Limits. 



Figure 74. Contents of Simulink block Mid Value. 
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Figure 75. Contents of Simuiink block Bias. 



ICvalue 


Figure 76. Contents of Simuiink block VarLimitlnteg. 


Figure 77 shows the contents of the block VTerrorWashout within the speed control law in 
Figure 71. The value for the gain KWO is two. Figure 78 shows the contents of block 
VertSpdControl within the speed control law . The limits for the limit block on the input signal 
to the integrator on the left side of Figure 78 are set to ±VSHddLim. 



Figure 77. Contents of Simulink block VTerrorWashout. 



VNAV_VS_Eng 


Figure 78. Contents of Simulink block VertSpdControl. 


Figure 79 shows the contents of block FadlnOut that smooths the pitch steering command signals 
when control mode switching of the path, speed, and vertical speed control laws occur and when 
switching to a new path segment occurs. 



Figure 79. Contents of Simulink block FadelnOut. 
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Figure 80 shows the contents of block PitchCmdProcess located at the upper right of Figure 47. 
The blocks in Figure 80 compute the pitch command (TF1ETAC) and the pitch error (DTHETE); 
then, DTHETE is used by the pitch inner loop control system that computes the elevator surface 
command. The pitch command is computed from the pitch steering command signals (VNAVS 
and VNAVI) that are outputs from the VNAV control system. Figure 81 shows the contents of 
block Subsystem. 



VNAVI 


Figure 80. Contents of Simulink block PitchCmdProcess. 



ICtrig 


Figure 81. Contents of Simulink block Subsystem. 
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